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ABSTRACT 

We present the results of a deep 300 ks Chandra HETG observation of the highly variable 
narrow-line Seyfert Type 1 galaxy NGC 405 1 . The HETG spectrum reveals 28 significant soft 
' X-ray ionised lines in either emission or absorption; primarily originating from H-like and He- 

like K-shell transitions of O, Ne, Mg and Si (including higher order lines and strong forbidden 
emission lines from O VII and Ne IX) plus high ionisation L-shell transitions from Fe XVII to 
Fe XXII and lower ionisation inner-shell lines (e.g. O VI). Modelling the data with XSTAR re- 
quires four distinct ionisation zones for the gas, all outflowing with velocities < 1 000 km s _1 . 
A selection of the strongest emission/absorption lines appear to be resolved with FWHM of 
^600 km s^ 1 . We also present the results from a quasi-simultaneous 350ks Suzaku observa- 
tion of NGC 405 1 where the XIS spectrum reveals strong evidence for blueshifted absorption 
lines at ^6.8 and ^7.1keV, consistent with previous findings. Modelling with XSTAR sug- 
CN| , gests that this is the signature of a highly ionised, high velocity outflow (log£ = 4.lio!i> 

v out ~ —0.02c) which potentially may have a significant effect on the host galaxy environ- 
ment via feedback. Finally, we also simultaneously model the broad-band 2008 XIS+HXD 
Suzaku data with archival Suzaku data from 2005 when the source was observed to have 
entered an extended period of low flux in an attempt to analyse the cause of the long-term 
spectral variability. We find that we can account for this by allowing for large variations in the 
O^l . normalisation of the intrinsic power-law component which may be interpreted as being due to 

significant changes in the covering fraction of a Compton-thick partial-coverer obscuring the 
central continuum emission. 
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1 INTRODUCTION outflowing winds (ranging in velocity from several hundred to sev- 
eral thousand kms , Blustin et al. 2005) and arise from a wide 

Recent systematic X-ray studies of AGN (Active Galactic Nuclei) range of column densities and levels of ionisation. Through the 

with ASCA, Chandra, XMM-Newton and Suzaku have established study of blueshifted absorption lines of K-shell Fe, recent observa- 

that at least half of all active galaxies host photo-ionised "warm" tions of higher luminosity AGN have also suggested the presence 

absorbers (e.g. Reynolds & Fabian 1995; Blustin et al. 2005) and of highly ionised ("hot") absorbers originating in high- velocity disc 

that they may in fact be ubiquitous in AGN, becoming observable winds (v out ~ —0.1c) with around a dozen cases to date (e.g. PG 

only in certain lines of sight (Krongold et al. 2008; although one 121 1+143, Pounds et al. 2003; PDS 456, Reeves, O'Brien & Ward 

other interpretation is that ~50 per cent of AGN do not have warm 2003) although there are a set of lower luminosity sources that 

absorbers). When observed at high spectral resolution (such as with still exhibit deep absorption lines at Fe K but require high column 

the transmission gratings on-board Chandra), these absorbers are densities and have outflow velocities on the order of a few thou- 

known to produce numerous narrow absorption features ranging sand kms -1 (e.g. Mrk 776, Miller et al. 2007; NGC 3516, Turner 

from various elements such as oxygen, silicon, neon, carbon, ni- et al. 2005; NGC 1365, Risaliti et al. 2005; NGC 3783, Reeves 

trogen, sulphur, magnesium and iron (e.g. Kaspi et al. 2002; Cren- et al. 2004). Recently, Tombesi et al. (2010) found evidence for 

shaw, Kraemer & George 2003; McKernan, Yaqoob & Reynolds blueshifted features in 17 objects within their sample of 42 radio- 
2007). The associated lines are often blueshifted, thus implying 
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quiet AGN observed with XMM-Newton, detecting 22 absorption 
lines at rest-frame energies >7.1 keV. This implies that high veloc- 
ity outflows may be a common feature in radio-quiet AGN. Such 
ultra-fast outflows indicate that the gas must be located very close 
to the nucleus of the AGN and so deep studies of these objects are 
vital to establish a greater understanding of the outflow kinematics 
and locations relative to the central continuum source which in turn 
can help to ultimately unravel the inner structure of AGN. 

One possible interpretation regarding accretion disc outflows 
is that they are perhaps produced as a result of radiation pressure 
a few R g from the event horizon (Proga, Stone & Kallman 2000). 
The high luminosities of these systems are the result of radiatively 
efficient accretion onto a supermassive black hole (SMBH) and 
such outflows could be a consequence of near-Eddington accre- 
tion (King & Pounds 2003). An alternative interpretation of many 
outflows is that they could be magnetohydrodynamically (MHD) 
driven as has been suggested for such objects as GRO J1655-40 
(Miller et al. 2008) and NGC 4151 (Kraemer et al. 2005; Crenshaw 
& Kraemer 2007). The derived outflow rates can in some cases be 
comparable to the mass accretion rate (several Solar masses per 
year) of the AGN and so in terms of kinetic power, the outflow 
can be responsible for a significant proportion of the bolometric lu- 
minosity. As a result, outflows are considered to be an important 
phenomenon in AGN and are believed to play a key role in feed- 
back processes between the black hole and the host galaxy (King 
2003, 2010), ultimately leading to the observed M-o relation for 
galaxies (Ferrarese & Merritt 2000; Gebhardt et al. 2000). 

NGC 4051 is a nearby (z = 0.002336; Verheijen & San- 
cisi 2001) narrow-line Type 1 Seyfert Galaxy at a distance of 
15.2 Mpc (Russel 2004) obtained from the Tully-Fisher relation for 
nearby galaxies (Tully & Fisher 1977). It has a black hole mass of 
M BH = l-73tg'.52 x 10 6 M Q determined via optical reverberation 
mapping (Denney et al. 2009) and is well known for its extreme 
X-ray variability both on short and long time-scales with the X-ray 
spectrum hardening as the source flux becomes lower (Lamer et al. 
2003), as seen in many other Type 1 Seyfert Galaxies (e.g. NGC 
5506, Lamer, Uttley & McHardy 2000; MCG-6-30-15, Vaughan & 
Edelson 2001). 

The warm absorber in NGC 4051 was studied by Collinge 
et al. (2001) with a simultaneous Chandra HETG (High Energy 
Transmission Grating) and Hubble Space Telescope (HST) Imag- 
ing Spectrograph observation where they discovered two separate 
X-ray absorption systems corresponding to outflow velocities of 
Wout = -(2340 ± 130) and v out = -(600 ± 130) kms" 1 . Nine 
separate absorption systems were detected in the HST UV spec- 
trum with one of the zones possibly corresponding to the lower 
velocity zone detected in the X-ray band. A further RGS (Reflec- 
tion Grating Spectrometer) observation with XMM-Newton was 
analysed by Ogle et al. (2004) who claimed the presence of a rel- 
ativistically broadened O VIII Lya emission line (EW ~ 90 eV) at 
~655 eV and suggested that this emission along with its associ- 
ated radiative recombination continuum (RRC) could be responsi- 
ble for the weak soft excess. However, they found no evidence of 
the high velocity outflow detected by Collinge et al. (2001). The 
RGS data from 2001 and 2002 were also analysed by Pounds et 
al. (2004a) and Nucita et al. (2010). They note that the 2002 ob- 
servation caught the source during an extended period of low flux 
CF2-10 = 5.8 x 10~ 12 ergcirT 2 s _1 ). Pounds et al. (2004a) con- 
cluded that the hard X-ray spectral shape observed during this pe- 



riod of low flux was due to an increase in opacity of a substantial 
column of gas in the line-of-sight causing the spectrum to become 
dominated by a quasi-constant cold reflection component. 

A more recent study by Krongold et al. (2007) showed the 
presence of two distinct ionisation components for the absorber 
through a lOOks XMM-Newton exposure. By measuring the elec- 
tron densities through the absorber variability, they inferred that 
the absorbing components must be compact with the high and low 
ionisation zones existing at distances of 0.5-1 1-d (light-days) and 
<3.5 1-d from the central engine respectively, well within the dusty 
torus and strongly suggestive of an accretion disc origin for the 
warm absorber wind. They calculated that the implied mass out- 
flow rate of the warm absorber wind corresponds to approximately 
2-5 per cent of the mass accretion rate of the source. They also 
detected several narrow emission lines in the RGS spectrum from 
C VI, N VI, O VII, O VIII, Ne IX and Fe XVII. 

A further study of the warm absorber was also undertaken 
by Steenbrugge et al. (2009) through high-resolution X-ray spec- 
troscopy using the Chandra LETG (Low Energy Transmission 
Grating). They discovered that this object contains an outflowing 
wind consisting of four separate absorbing components ranging in 
ionisation parameter, log£, from 0.07 to 3.19. They found that the 
absorbing zone with log £ = 3.19 requires a high outflow velocity 
of v out = — 4760 km s -1 and from a study of the variability of the 
warm absorber, they inferred that three of the four absorbing zones 
appear to be located in the range of 0.02-1 pc from the black hole. 

NGC 4051 was observed in 2005 with Suzaku when it was 
found to have fallen into an extended period of historically low 
flux (Terashima et al. 2009). A strong excess of emission was 
seen at energies >10keV suggesting that the primary power-law 
continuum had largely disappeared and that the bulk of the ob- 
served emission was reflection-dominated (e.g. Guainazzi et al. 
1998; Pounds et al. 2004a). However, Terashima et al. (2009) were 
able to describe the data using a partial covering model whereby the 
reflection-dominated emission and the intrinsic power-law emis- 
sion are independently absorbed by gas covering some significant 
fraction of 47r sr. In this scenario, changes in the covering fraction 
of the partially-covered power-law component can account for the 
changes in the spectral shape at low energies whereas changes in 
the normalisation of the power-law component overlaid on a nearly 
constant hard component can account for the variability at energies 
>3.5keV. This interpretation of the long-term spectral variability 
being caused by the primary X-ray continuum disappearing from 
view leaving behind a constant, hard, reflection-dominated compo- 
nent was supported by the results of Miller et al. (2010) through 
principal components analysis (PCA) of the 2005 and 2008 Suzaku 
data. Miller et al. (2010) also suggested that the low-flux states are 
caused by variable partial covering obscuring the central engine 
which also explains the constancy of the narrow FeKa emission 
line flux. 

Finally, the presence of an additional absorption feature at 
^7.1 keV was reported by Pounds et al. (2004a) indicating that 
a highly ionised, high velocity outflow could be apparent in this 
source. Further observations with Suzaku and XMM-Newton have 
confirmed the presence of this feature along with an additional 
absorption line at ~6.8keV (Terashima et al. 2009; Pounds & 
Vaughan 2010, submitted). If the two absorption lines are associ- 
ated with K-shell transitions from Fe XXV and Fe XXVI then their 
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blueshift would imply that the material is outflowing with a veloc- 
ity of w out ~ —0.02c. 

Here we report on a ~300ks Chandra HETG (High Energy 
Transmission Grating) observation of NGC 4051 where we aim 
to study the soft X-ray warm absorber in detail (Sections 4 and 
5). We also report on a contemporaneous 350 ks Suzaku observa- 
tion where we study the Fe K band including the blueshifted Fe K 
absorption lines as the signature of the highest ionisation compo- 
nent of the outflowing wind. We then proceed to analyse the en- 
tire broad-band Suzaku spectrum from 0.5-50.0 keV. Finally, with 
a combined analysis of the 2005 archival Suzaku data (exposure 
time ~ 120 ks) when the source was found to be in an extended 
period of low flux, we also aim to study the origin of the long-term 
X-ray variability of this AGN (Section 6). 



o 

o 




2 OBSERVATIONS AND DATA REDUCTION 
2.1 Chandra Analysis 

NGC 4051 was observed with the Chandra X-ray Observatory 
(Weisskopf et al. 2000) on twelve different occasions dating from 
2008 Nov 6 to 2008 Nov 30 with a total exposure time of ~300ks. 
A log of the separate observations and their corresponding expo- 
sures is shown in Table 1. The Chandra observations were per- 
formed with the High Energy Transmission Grating Spectrome- 
ter (HETGS; Markert et al. 1994; Canizares et al. 2005) in the 
focal plane of the Advanced CCD Imaging Spectrometer (ACIS; 
Garmire et al. 2003) which consists of two separate gratings: the 
Medium Energy Grating (MEG; 0.4-8.0 keV) and the High Energy 
Grating (HEG; 0.7-10.0 keV). The Chandra HETG data were re- 
duced using version 4.1 of both the Chandra Interactive Analysis 
of Observations software package (CIAO 1 ; Fruscione et al. 2006) 
and corresponding Calibration Database (CALDB). 

2.1.1 HETG Reduction 

First order MEG and HEG spectra were extracted for the source 
and background for each individual observation. Spectral redis- 
tribution matrix (RMF) files were created using the MKGRMF 
script for each first order grating arm (-1 and +1) for the MEG 
and HEG. Telescope effective area files were also created us- 
ing the FULLGARF script incorporating the CIAO tool MKGARF. 
Events files were extracted from the negative and positive first or- 
der grating arms of the MEG and HEG allowing light curves and 
spectra to be extracted. The first order spectra were then com- 
bined for each individual exposure using combined response files 
(with appropriate weighting) for the MEG and HEG. The back- 
ground was not subtracted however as it has a negligible effect 
in the energy ranges of interest here. The total count rates in 
the first-order energy spectra (0.5-10.0 keV) are 0.602 ± 0.001 
and 0.250 ± 0.001 cts~ x for the MEG and HEG respectively (or 
0.583 ± 0.001 and 0.233 ± 0.001 ct s" 1 for the respective 0.5-5.0 
and 1.0-8.0 keV energy bands for the MEG and HEG which we 
adopt here). These count rates correspond to fluxes of 3.99 x 10~ n 
and 3.29 x 10 _11 erg cm -2 s _1 (see Table 2) and luminosities of 
4.78 x 10 41 and 3.95 x 10 41 ergs _1 for the respective 0.5-5.0 
and 1.0-8.0keV energy bands for the MEG and HEG. A lightcurve 
showing the twelve individual observations with 2 ks binning over 

1 http://cxc.harvard.edu/ciao 



Figure 1. The summed HEG and MEG lightcurve showing the count rates 
of the twelve individual Chandra HETG observations of NGC 405 1 over the 
1.0-6.0 keV energy range. The labels refer to the sequence numbers of the 
observations listed in Table 1 . The data are binned up in 2 ks bins and the 
horizontal lines show the duration and overlap of the two contemporaneous 
Suzaku observations. 

the 1 .0-6.0 ke V range (including both the MEG and HEG) is shown 
in Figure 1 . 

2.2 Suzaku Analysis 

NGC 4051 was observed with Suzaku in November 2005 and 
again in November 2008 with net exposures of 120 and 350 ks 
respectively 2 (see Table 1 for an observation log). The 2005 data 
have been described previously by Terashima et al. (2009). Here 
we discuss data taken with the Suzaku XIS (X-ray Imaging Spec- 
trometer; Koyama et al. 2007), consisting of four X-ray telescopes 
(Mitsuda et al. 2007) each with a CCD in the focal plane, and the 
PIN diodes of the non-imaging HXD (Hard X-ray Detector; Taka- 
hashi et al. 2007). Note that NGC 4051 is too faint to be detected 
with the HXD GSO instrument. Events files from versions 1.2 and 
2.2.1 1.22 of the Suzaku pipeline processing were used for the 2005 
and 2008 observations respectively. All data were reduced using 
HEASOFT version 6.4.1. 

2.2.1 XIS Reduction 

A full description of the Suzaku XIS data reduction procedure can 
be found in Turner et al. (2010) since the same base files are used 
here. The count rates corresponding to the 2005 Nov 10, 2008 Nov 
6 and 2008 Nov 23 observations respectively over the 0.5-10.0 keV 
energy range were found to be ~0.45, ~2.07 and ~1.42cts _1 per 
front-illuminated XIS with the background rates corresponding to 
~2.4 per cent, ~0.9 per cent and ~0.9 per cent of the respec- 
tive source count rates. These count rates correspond to fluxes of 
1.32 x 10" 11 , 5.03 x 10~ n and 3.53 x 10~ n erg cm -2 s" 1 and 
luminosites of 3.66 x 10 41 , 1.57 x 10 42 and 1.10 x lO^ergs -1 
for the three respective observations. The three front-illuminated 
XIS chips (XIS 0, 2, 3; hereafter XIS-FI) were used in the spectral 

2 Note that the 2008 Suzaku observation was split into two separate obser- 
vations due to scheduling constraints with respective exposure times of 275 
and 78 ks. 
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Table 1. Suzaku and Chandra HETG observation log of NGC 4051. 



Mission Sequence Number 



Observation Date 
(Start Date) 



Exposure (ks) 



Suzaku 


700004010 


2005-11-10 


120 


Suzaku 


703023010 


2008-11-06 


275 


Chandra 


10777 


2008-11-06 


27.8 


Chandra 


10775 


2008-11-08 


30.9 


Chandra 


10403 


2008-11-09 


38.2 


Chandra 


10778 


2008-11-11 


34.1 


Chandra 


10776 


2008-11-11 


25.2 


Chandra 


10404 


2008-11-12 


20.1 


Chandra 


10801 


2008-11-13 


26.1 


Chandra 


10779 


2008-11-20 


27.8 


Suzaku 


703023020 


2008-11-23 


78 


Chandra 


10780 


2008-11-25 


26.4 


Chandra 


10781 


2008-11-26 


24.1 


Chandra 


10782 


2008-11-29 


23.7 


Chandra 


10824 


2008-11-30 


9.2 



We also used the cross-sections for X-ray absorption by the inter- 
stellar medium obtained by Morrison & McCammon (1983) by us- 
ing WABS in XSPEC. Abundances are those of Anders & Grevesse 
(1989) unless otherwise stated. Due to the low number of counts per 
channel (i.e. <50 counts per resolution bin) in the HETG data, we 
could not use \ 2 minimisation; instead, we checked our fit statis- 
tics in the Chandra spectra using the C-statistic (Cash 1979). How- 
ever, due to the higher photon statistics in the Suzaku spectra, all 
goodness-of-fit values were checked using \ 2 minimisation. Note 
that all errors quoted correspond to 90 per cent confidence for one 
interesting parameter (AC = Ax 2 = 2.71) unless stated other- 
wise. Where AC or A\ 2 values are quoted in tables, these values 
have been determined by removing the component from the final 
model and re-fitting. All fit parameters are given in the rest frame of 
the host galaxy having been corrected for the cosmological redshift 
(z = 0.002336) and assuming a distance of 15.2 Mpc to NGC 405 1 
(Russell 2004) obtained from the Tully-Fisher relation for nearby 
galaxies (Tully & Fisher 1977). 



analysis of the 2005 data due to their greater sensitivity at Fe K. 
However, as the use of XIS 2 was discontinued after a charge leak 
was discovered in 2006, the spectral analysis from the 2008 data 
refers only to XIS and 3. The XIS 1 (back-illuminated) data 
were checked for consistency but not used during the spectral fit- 
ting due to the detector's reduced sensitivity compared to XIS-FI 
at FeK. As the chips were found to produce consistent spectra 
within the statistical errors in all observations, the XIS-FI spec- 
tra and responses were combined to maximise signal to noise. The 
XIS source spectra were binned at the half-width at half-maximum 
(HWHM) resolution of the detector due to the high photon statis- 
tics. This enabled the use of \ 2 minimisation in all Suzaku fits as 
there were >50 counts per resolution bin. 

2.2.2 HXD Reduction 

Likewise, the Suzaku HXD data reduction procedure can also be 
found in Turner et al. (2010). As NGC 4051 is too faint to be de- 
tected by the HXD GSO instrument, we use data taken with the 
HXD PIN only which provides useful data over the 15.0-70.0 keV 
energy range. Our net exposure times of the PIN source spectra 
were found to be 1 12, 204 and 59 ks for the 2005 Nov 10, 2008 Nov 
6 and 2008 Nov 23 observations respectively. Note that the HXD 
PIN response file dated 2008/01/29 (epoch 1) for the HXD nom- 
inal pointing position was used for the 2005 observation whereas 
the response file dated 2008/07/16 (epoch 5) for the XIS nominal 
pointing position was used for the the 2008 observations. For the 
three respective observations, the net PIN source count rates are 
0.043 ± 0.001, 0.068 ± 0.001 and 0.053 ± 0.001 cts" 1 over the 
15-50 keV range. These count rates in turn correspond to fluxes of 
1.73 x 10 -11 , 3.06 x 10 -11 and 2.66 x 10 -11 erg cm" 2 s" 1 for the 
three respective observations. The mean net count rates and fluxes 
are summarised in Table 2. 



3 SPECTRAL ANALYSIS 

For all spectral analysis we used the XSPEC V 1 1 .3 software pack- 
age (Arnaud 1996). All our fits include Galactic absorption with a 
column density N§ AL = 1.35 x 10 20 cm~ 2 , obtained using the 
FTOOL NH with the compilations of Dickey & Lockman (1990). 



4 THE CHANDRA HETG SPECTRUM AT LOW 
ENERGIES 

We initially considered the twelve separate HETG observations 
(see Table 1) and found that the flux appears to vary by a factor of 
~2 over the course of this 24-day period. The short-term variabil- 
ity appears to be consistent with bright narrow-line Type 1 Seyfert 
Galaxies in that the X-ray continuum softens as the flux increases 
(e.g. Lamer et al. 2003; Taylor, Uttley & McHardy 2003). How- 
ever, we find that the low signal-to-noise ratio due to these short 
individual exposures is insufficient to significantly determine the 
parameters of the warm absorber and detect individual absorption 
lines so we begin by considering the entire 300 ks time-averaged 
HETG spectrum. We initially binned the MEG and HEG data to a 
constant resolution of AA = 0.01 A and AA = 0.005 A respec- 
tively and in all subsequent fits we consider the MEG and HEG 
over the respective 0.5-5.0 keV and 1 .0-8.0 keV energy bands. 

4.1 Continuum 

We initially fitted a simple parameterisation of the continuum us- 
ing the MEG and HEG in order to allow the warm absorber to be 
studied in detail. The continuum emission was fitted over the 0.5- 
8.0 keV range with an absorbed power-law with a Galactic column 
density N§ AL = 1.35 x 10 20 cm" 2 (Dickey & Lockman 1990) 
and a photon index T = 2.12 ± 0.01. This gave a poor fit to the 
data with C/d.o.f. = 6082/4396. The fluxed HETG spectrum un- 
folded against a power-law with V = 2 is shown in Figure 2 where 
significant emission can be observed >6 keV due to the Fe K emis- 
sion complex and positive spectral curvature can also be seen below 
1 keV indicating that an additional steepening of the spectrum is re- 
quired in this source. We parameterised this with the addition of a 
simple featureless black body (kT — 0.10 ± 0.01 keV) which im- 
proved the fit by AC = 197. The purpose of the black body com- 
ponent was not to model the data in a physical sense but to simply 
parameterise the soft excess with a smooth continuum, therefore al- 
lowing the individual lines from the warm absorber to be studied in 
detail. This yields a decent fit to the low-energy continuum and as 
we see no evidence for any emission / absorption features >2keV 
(apart from at Fe K), we excluded the MEG and HEG data >2 keV 
and proceeded to analyse and identify the soft ionised absorption 
lines in the 0.5-2.0 keV band. A plot of the a residuals <2keV is 
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Table 2. Mean count rates and fluxes for the time-averaged Chandra HETG and Suzaku observations of NGC 4051. All XIS, HXD, MEG and HEG count 
rates and fluxes are given in the 0.5-10.0, 15.0-50.0, 0.5-5.0 and 1.0— 8.0 keV energy bands respectively. 



Mission 



XIS 

Suzaku m „ 2005-11-10 



HXD 
XIS 
HXD 
MEG 
HEG 
XIS 
HXD 



XTS 

Suzaku J™ 2008-11-06 



Cnandra ***** 2008-11-06 



XTS 

Suzaku 2008-11-23 



Count Rate 


Flux 


ICl S J 


^A_LU Clg Clll 


0.45 


1.32 


0.04 


1.73 


2.07 


5.03 


0.07 


3.06 


0.58 


3.99 


0.23 


3.29 


1.42 


3.53 


0.05 


2.66 



- 1 ) 




0.5 1 2 5 



Observed Energy (keV) 

Figure 2. A plot of the HETG spectrum from 0.5-8.0 keV unfolded against 
a power-law with Y = 2. Strong emission at energies >6 keV is most likely 
due to significant FeK emission and positive spectral curvature can also 
be observed at energies < 1 keV suggestive of the presence of a weak soft 
excess. The MEG and HEG are fitted over the 0.5-5.0 and 1.0-8.0 keV 
bands respectively (in the observed frame) and are binned by a factor of 10 
for clarity. 

shown in Figure 3 revealing a number of absorption and emission 
lines. 



4.2 Absorption Lines 

The time-averaged 300 ks spectrum reveals 21 significant (>99.9 
per cent) soft X-ray absorption lines below 2.0 keV originating 
from material with a range of ionisation states and column den- 
sities. Several of the lines appear to be somewhat resolved in the 
HETG spectrum (e.g. O VII ls-2p, O VIII ls-2p, Fe XVIII 2p-3d) 
with FWHM on the order of a few hundred to ~1 000 km s -1 . In- 
deed we find a mean value for the FWHM of 580 ± 59 kms - 
with a dispersion of a = 210 km s -1 for the 13 absorption lines 
which appear to be resolved in the spectrum. The spectrum re- 
veals absorption primarily from H-like and He-like ions of oxygen, 
neon, magnesium and silicon along with lines from less-ionised 
ions such as Ovi. We also detect several L-shell transitions of 
Fe XVII-Fe XXII. After taking into account the systemic redshift of 
NGC 4051, we note that the observed centroid energies of the ab- 
sorption lines are all blueshifted on the order of a couple of eV. 
This blueshift implies that the absorbing material in our line-of- 



sight is outflowing with calculated velocity shifts on the order of a 
few hundred to ~1 000 km s -1 (e.g. the O VII, O VIII ls-2p lines) 
relative to the host galaxy. A mean value for the outflow velocity of 
Vout = —(620 ± 34) km s _1 with a dispersion of a = 150 kms -1 
is calculated for the 20 absorption lines that we are able to identify. 
We determined the parameters of the lines by modelling them with 
simple Gaussian profiles. These parameters and their likely iden- 
tifications are summarised in Table Al. All lines were detected at 
the >99.9 per cent significance level for two interesting parameters 
(i.e. AC > 13.8). 



4.3 Emission Lines 

In addition to a wealth of absorption features, the HETG spectrum 
also reveals 7 narrow emission lines which were again initially 
modelled with simple symmetric Gaussians. Their identifications 
appear to be largely consistent with the narrow emission lines found 
in the XMM-Newton RGS spectrum of NGC 4051 by Pounds et al. 
(2004a) and also with those detected by Terashima et al. (2009). 
These lines appear to primarily originate from He-like ions of oxy- 
gen, neon and silicon due to their associated forbidden transitions. 
The NelX intercombination lines are also identified at ~915 eV 3 . 
Furthermore, Lya transitions from H-like species of oxygen and 
magnesium are also observed at ^653 eV and ~1472eV respec- 
tively. We note that the O VIII Lya emission that we observe here 
appears to be narrow (FWHM <400kms _1 ; see Figure 4) and 
that we do not require an additional relativistically broadened Lya 
feature such as the one claimed by Ogle et al. (2004). The emis- 
sion lines that we detect here have a mean velocity shift, v out , 
within ±200 km s -1 of the systemic shift and appear to be unre- 
solved except for the forbidden emission lines from Ne IX (FWHM 
= 550t 3 ™ kms" 1 ) and Si XIII (FWHM = 1 000± 7 4 9 ° km s" 1 ). We 
do note that the Lya transitions of O VIII and Mg XII appear to be 
redshifted on the order of a few 100 kms -1 relative to their rest 
frame centroid energies. This is most likely due to the P-Cygni-like 
profiles in which the lines appear with their corresponding absorp- 
tion lines (see Figure 4). The best-fitting parameters of the emission 
lines are listed in Table Al. Like the absorption features, all of the 
emission lines were also detected at the >99.9 per cent significance 
level for two interesting parameters (i.e. AC > 13.8). 



3 Interestingly, the forbidden and intercombination lines from Ne IX appear 
to have similar strengths; this could be due to the presence of the Ne VIII 
ls-2p absorption line at ~895 eV reducing the observed flux of the Ne IX 
forbidden transition, thus making it appear to be weaker than it really is. 
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Figure 3. The <r residuals of the HETG spectram from 0.5-2.0 keV when fitted with a single absorbed power-law and a featureless blackbody as described in 
Section 3.1. Significant residuals can be observed throughout the spectrum indicating the presence of several absorption and emission features. The statistically 
significant features which we are able to identify are labelled on the plot. All spectra are binned up by a factor of 2 for clarity and are plotted in the observed 
frame. 
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5 THE BROAD-BAND HETG SPECTRUM 

In order to model the broad-band HETG spectrum we included the 
data above 2keV in our spectral fits so that we were considering 
the full 0.5-8.0 keV energy range. We carried forward the model 
from Section 3 consisting of a simple baseline continuum, a fea- 
tureless blackbody to parameterise the soft spectral steepening and 
symmetric Gaussians to model the soft absorption and emission 
features. Figure 2 shows that a significant excess can be observed 
in the FeK band at energies >6keV (see Section 5.2). However, 
we begin by considering a more physical parameterisation of the 
warm absorber. 



5.1 The Warm Absorber 

We removed the simple Gaussian profiles parameterising the ab- 
sorption lines from our model and instead attempted to model 
the warm absorber using the XSTAR 2.11nll code of Kallman & 
Bautista (2001; also see Kallman et al. 2004), incorporating the 
abundances of Grevesse, Noels & Sauval (1996). The XSTAR code 
assumes thin shells of absorbing gas and self-consistently models 
zones of absorption parameterised by their column density, Nb_, 
and ionisation parameter, f , which is defined as: 



(1) 



and has units erg cms - , where Li on is the ionising luminos- 
ity from 1 to 1 000 Rydbergs in units erg s , n is the gas density 
in cm -3 and R is the radial distance in cm of the absorbing gas 
from the central source of X-rays. For the spectral energy distribu- 
tion (SED) of the XSTAR models we assume a simple illuminating 
power law with T — 2.5, although this is difficult to compare with 
the observed SED of NGC 4051 since the UV data from the HST 
were acquired a year later and so are not simultaneous. We initially 
attempted to model the absorption lines using Solar abundances and 
a turbulent velocity width of a = 200 km s ~~ 1 , a value largely com- 
parable to the observed velocity widths from the Gaussian line fits 
(see Table Al). 

We find that we statistically require four individual zones of 
absorbing gas to model the data covering four distinct levels of 
ionisation. Each zone is significant at the >99.9 per cent confi- 
dence level. The column densities are found to be on the order of 
JVh ~ 10 20 -10 21 cm -2 and the zones appear to cover a range in 
ionisation parameter of log £ = —0.86 to log£ = 2.97. All the 
zones appear to be blueshifted implying that the material is out- 
flowing with velocities on the order of a few 100km s _1 with the 
general trend being that v ou t increases with increasing £ . We find 
that one of the zones requires two significantly different outflow 
velocities and so an additional zone with comparable values for 
the ionisation parameter and column density is included in the fit 
to account for this (see zones 3a and 3b in Table 3). We note that 
replacing zones 3a and 3b with one zone of higher turbulence ve- 
locity significantly worsens the fit. Therefore, the inclusion of the 
two separate zones with a = 200 km s -1 and differing outflow 
velocities appears to be preferred by the data. Note that we subse- 
quently refer to zones 3a and 3b as one single zone. A summary of 
the best-fitting values and the corresponding fit statistic (AC) for 
each zone is given in Table 3. 
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Figure 4. A plot showing the relatively narrow O VIII Lyo emission at 
~652eV in the observed frame. This line, along with its corresponding 
absorption feature at ~653 eV appears to be P-Cygni-like in appearance. 
The model is superimposed on the data and is shown in red; the O VIII hya 
emission at ~652 eV is modelled with a Gaussian and the corresponding 
absorption line at ~653 eV is modelled with XSTAR as described in Section 
5.1. We note that a turbulence velocity of a = 500 km s — 1 in the warm ab- 
sorber model (Zone 2) is required to model the absorption line at ~666 eV 
which likely corresponds to the ls-3p transition from O VII. 



Regarding the transitions modelled by each XSTAR zone, we 
find that the first zone (ultra-low ionisation) is mainly responsible 
for the absorption of O VI and lower ionisation O ions. The sec- 
ond zone, with a slightly higher ionisation parameter, appears to 
model absorption primarily from O VII. An inspection of the resid- 
uals reveals that this zone under-predicts the amount of absorp- 
tion at ~666 eV which is likely due to the turbulence velocity of 
a — 200 km s" 1 being too low. This has the effect of saturating the 
O VII ls-3p line and hence underpredicting its EW. Increasing the 
turbulence velocity of this zone to a — 500 km s -1 significantly 
improves the fit and better models the residuals at ~666eV (see 
Figure 4). However, we find that increasing the turbulence velocity 
of the remaining warm absorber zones to a — 500 km s -1 signifi- 
cantly worsens the fit and that these zones are better modelled with 
a — 200 kms -1 . Zones 3a and 3b appear to adequately model the 
absorption from ions such as Ne IX, O VIII and the various L-shell 
transitions from Fe XVII - XXII and finally, the highest ionisation 
zone appears to correspond to absorption from ions such as Ne X, 
Mg XI - XII and Si XIII. A plot of the warm absorber model super- 
imposed on the data is shown in Figure 5 and a further plot show- 
ing the individual contribution of each of the four zones is shown 
in Figure 6. A subsequent paper (Crenshaw et al, 2010 in prepa- 
ration) will discuss the results of a UV observation of NGC 4051 
taken with the Cosmic Origins Spectrograph (COS) on-board the 
HST. The preliminary results suggest that the X-ray zones 1 and 2 
detected here with Chandra may also coincide with absorption in 
the UV data. 

We note there are no detected absorption lines associated 
with a fifth warm absorber zone. However, we did allow for the 
presence of a partially-covering zone of absorption by including 
an additional power-law component absorbed by a further shell 
of gas. We tied the photon index of the two power-law compo- 
nents together leaving just one of the components absorbed by 
an XSTAR zone. We note that zones 1-4 of the warm absorber 
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and a column of Galactic hydrogen remained in the model ab- 
sorbing the entire continuum. We find that this additional partial- 
coverer improves the fit significantly with AC = 286 and appears 
to be important in modelling additional curvature in the contin- 
uum at energies below about 5 keV. The best-fit parameters of the 
partial-coverer indicate that the material is partially-ionised with 
log£ < 2.40 and has a substantial line-of-sight column density, 
N H = 2.00±S:|1 x 10 23 cm" 2 . The addition of this component 
also allows the temperature of the black body component to re- 
adjust slightly to kT = 0.12 ± 0.01 keV (a value consistent with 
Terashima et al. 2009) and the best-fitting photon index of the 
power-law continuum to steepen to V = 2.25+° ) ' ) i- Indeed, with- 
out the partial-coverer, we note that the photon index of T ~ 2.1 
appears to be slightly too flat for this source in the high-flux state 
(e.g. Lamer et al. 2003; Miller et al. 2010). The normalisation 
ratio of the absorbed to total power-law components at 1 keV is 
4.98 x 10~ 3 /1.67 x 10" 2 w 0.30 which suggests that the partial- 
coverer may correspond to a line-of-sight covering fraction of ~30 
per cent of the nuclear X-ray source. The best-fitting parameters of 
the partial-coverer are also noted in Table 3. 

Furthermore, replacing the black body with the REFLIONX 
reflection model of Ross & Fabian (2005) provides an equally 
good parameterisation of the soft excess. We find that no rela- 
tivistic blurring is required and that the addition of this compo- 
nent allows the photon index of the power-law to steepen further 
with r = 2.45 ± 0.02, a value consistent with that found from 
a quasi-simultaneous broad-band Suzaku observation in 2008 (see 
Section 6 where this is discussed in more detail). The REFLIONX 
component has a best-fitting value for the ionisation parameter of 
log ^ = 3.3 ±0.1 and corresponds to a ratio of reflected flux to 
incident flux of ~0.8 over the 0.5-8.0 keV energy range. However, 
the ionised reflector is unable to account for some of the soft X-ray 
emission lines, particularly the forbidden transitions from Ovil, 
Ne IX and Si XIII which indicates that these lines could have an al- 
ternative origin. There have been previous claims in the literature of 
detections of RRC associated with H-like and He-like ions of C, N, 
O and Ne (e.g. Ogle et al. 2004; Pounds et al. 2004a) using the RGS 
on-board XMM-Newton which suggests that the soft X-ray emis- 
sion lines could instead be the signature of a photo-ionised plasma. 
Indeed it has been somewhat well established that such soft X-ray 
lines are commonly formed either in the narrow-line region (NLR) 
or as part of the photo-ionised outflow (e.g. Bianchi, Guainazzi & 
Chiaberge 2006). 



5. 7. 7 Comparison with Recent RGS Data 

The warm absorber in NGC 4051 was also recently studied by 
Pounds & Vaughan (2010) through XMM-Newton RGS observa- 
tions in May-June 2009. The data revealed significant absorption 
features primarily from H-like and He-like ions of C, N, O, Ne, 
Mg, Si, Ar and Fe. When modelled with XSTAR, the RGS data ap- 
pear to require five separate line-of-sight outflow velocity compo- 
nents ranging in velocity from v out ~ — 500 km s -1 to v out ~ 
-30 000 km s _1 . With our HETG data, we are able to confirm here 
the presence of the lowest velocity zone (on the order of a few 
hundred kms -1 ) and also the third zone (v out ~ -6 000 km s" 1 ) 
which can clearly be seen at Fe K with the Suzaku XIS (see Section 
6.3 for more details). However, we find no requirement for the high- 
est velocity zone at FeK (« ou t ~ — 30 000 km s -1 ; Section 6.3); 
nor do we find any evidence for the second and fourth intermediate 
zones claimed by Pounds & Vaughan (2010) which have velocities 
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Figure 6. Plot showing the contribution of each of the four warm absorber 
zones to the power-law continuum. The lowest to highest ionisation zones 
(zones 1 to 4) are shown from the upper to lower panels respectively. The 
ions contributing to the most significant absorption lines are marked on the 
plot. 



Table 3. The best-fitting parameters for the four individual zones of absorp- 
tion in the HETG data. a XSTAR zone: Nn, column density; ?, ionisation 
parameter; v ou t , outflow velocity. b Partial-coverer: / CO v , covering fraction. 
Zones 1, 3 and 4 require a turbulent velocity width of <j = 200 km s -1 
whereas Zone 2 is better modelled with a turbulent velocity width of 
<j = 500 kms" 1 . 
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respectively. 



-4 000 km s" 1 and v Q 
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Pounds & Vaughan (2010) report that the Ovil ls-2p reso- 
nance transition appears to be detected with three separate veloc- 
ity components (-400kms _1 ; -4000kms _1 ; -6 000 km s" 1 ). 
However, we are only able to confirm the lowest velocity compo- 
nent with the HETG for this transition (see Figure 7; top panel). 
Pounds & Vaughan (2010) also report four separate velocity com- 
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Figure 5. The HETG spectrum from 0.5-2.0 keV with our best-fitting warm absorber model super-imposed. The emission lines are modelled with simple 
Gaussians. All spectra are binned up by a factor of 2 for clarity and are plotted in the observed frame. Note that an additional Gaussian with a negative flux 
was included in the fit to model the excess residuals at ~666 eV (see Section 5.1 for further details). We do note that some unmodelled residuals remain in the 
data at ~0.87 keV; these are likely due to inner-shell Ne transitions (e.g. Ne VIII) which are not included in the XSTAR models. 
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Figure 8. Ratio plot showing the residuals of the HEG data from 5.5— 
7.5 keV to a power-law absorbed by Galactic hydrogen and a partial-coverer 
(as described in Section 5.1). Significant FeK line emission is clearly 
present. The data are binned by a factor of 2 for clarity. The vertical dot- 
ted lines show the expected line energies of, from left to right, Fe I-XVII Ka 
(red), Fe XXV forbidden (blue), Fe XXV Ka resonance (cyan), Fe XXVI Lya 
(magenta) and Fe I-XVII K/3 (orange) in the observed frame. 



Figure 7. Upper panel: Plot showing a radial velocity profile centred on the 
O VII ls-2p resonance transition at ~574 eV. Note that we see evidence 
for an associated absorption component with an outflow velocity of a few 
hundred kins -1 but we find no requirement for the higher velocity com- 
ponents reported by Pounds & Vaughan (2010). Lower panel: Plot showing 
a radial velocity profile centred on the O VIII Lya transition at ^653 eV. 
Again, we see evidence for the lowest velocity absorption component with 
«out < l,000kms -1 .Wedo see an absorption feature which may corre- 
spond to O VIII Lya blueshifted by ti ou t ~ —6 000 km s -1 ; however, we 
may also be able to associate this feature with the ls-3p transition from 
OVII. 



ponents for the O VIII Lya absorption line. However, we are again 
only able to confidently confirm the lowest velocity component in 
the HETG data (see Figure 7; lower panel). Interestingly, we do ob- 
serve additional residuals at ~666 eV in the HETG spectrum which 
could potentially be the signature of a higher velocity component 
of the O VIII Lya transition with w ou t ~ —6 000 km s _1 (Figure 7; 
lower panel), consistent with the value of v out ~ — 5 600 km s - 
found in the RGS spectrum (we note that we can confirm this veloc- 
ity component at Fe K with the Suzaku XIS; see Section 6.3). How- 
ever, as this residual feature coincides with the expected energy of 
the ls-3p transition from O VII, we may also be able to account for 
this with a slightly higher turbulence velocity (a = 500 km s _1 ) in 
our XSTAR model (as mentioned in Section 5. 1) so its interpretation 
remains open. We also note that we find no requirement at all here 
for the further two velocity components associated with the O VIII 
transition at v out ~ —4 100 km s _1 and v out ~ —9 000 km s _1 in 
the RGS data. 

Furthermore, Pounds & Vaughan (2010) also report that the 
Ne X and Mg XII Lya absorption lines are blueshifted on the order 
of f ut ~ — 7 000 km s -1 in the RGS data whereas we only find 
evidence for these components with v ou t < 1 000 km s _1 with the 
Chandra HETG. Therefore, we are only able to confirm two of the 
five velocity components associated with the XMM-Newton data, 



although we cannot rule out the possibility that the absorber may 
have significantly varied between observations. 

5.2 The FeK Complex 

A ratio plot of the residuals to the absorbed baseline HETG con- 
tinuum from 5.5-7.5 keV is shown in Figure 8. Line emission is 
clearly present suggestive of a significant Fe K complex. The most 
prominent line appears to be the FeKa fluorescence line from 
near-neutral material at En nc — 6.41 ± 0.03 keV. Modelling this 
with a Gaussian reveals that this line has an intrinsic width of 
a = 0.12±g ^keV (FWHM = 13 000^ 000 kms" 1 ) and an as- 
sociated equivalent width of EW = 1851 ?J eV. The addition of 
this line significantly improved the fit by AC = 99. However, 
the fit was improved further by incorporating an additional un- 
resolved narrow component (a = OeV) centered at 6.40 keV to 
model any Fe Ka emission from distant material. This component 
has an equivalent width of EW = 50j^g eV and resulted in an 
improvement to the fit statistic of AC = 19. Upon the addition 
of this component, the parameters of the underlying broader Ka 
component adjusted slightly to compensate but were still consis- 
tent within the errors (see Table 4 for the final best-fitting values). 
We did also test for the presence of a Compton shoulder by includ- 
ing an additional Gaussian centred around ~6.3 keV. This feature is 
statistically unrequired by the data and an upper limit on the equiv- 
alent width is found to be EW < 6 eV. Including this component 
did not result in any changes to the parameters of the broad and nar- 
row Fe Ka components at ~6.4 keV. So although the near-neutral 
Fe Ka emission appeared to then be well modelled, further resid- 
uals were still present indicating additional K-shell emission from 
ionised Fe. 

We adopted additional Gaussians to model the remaining 
emission. A line at En na = 6.97 ± 0.02 keV is required by 
the data with an intrinsic width of a = 32±f* eV (FWHM = 
3 100 "" 



63 



+32 



■ l 600 ' fln s X ) w< ^ an associated equivalent width of EW = 
eV. This line is most likely due to the ls-2p doublet from 
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hydrogen-like iron (Fe XXVI) and corresponds to an improvement 
in the fit statistic of AC = 15. A third component with a centroid 
energy of -Eii nc = 6.64±0.02 keV is also required by the data. This 
line has an intrinsic width of a < 52 eV (FWHM < 5 400 km s _1 ) 
and corresponds to an equivalent width of EW = 18 j] J| eV. This 
line improves the fit statistic by AC = 6 (i.e. only ~95 per cent 
significant) and is most likely associated with being the forbidden 
transition of the helium-like iron (Fe XXV) triplet at 6.636 keV, sug- 
gesting that this emission line may be consistent with an origin in 
photo-ionised gas (Bautista & Kallman 2000; Porquet & Dubau 
2000). No further emission or absorption features were observed 
around the FeK complex in the HETG data and the best-fitting 
values of the emission lines are shown in Table 4. The overall fit 
statistic corresponds to C/d.o.f. = 4779/4426. 

With the parameterisation of the baseline continuum, the 
warm absorber and the FeK components complete, we consider 
this to be our final best-fitting model to the time-averaged Chan- 
dra HETG data. We then proceed to apply this model to the 2008 
contemporaneous Suzaku data. 



6 SUZAKU SPECTRAL ANALYSIS 

We begin by considering the time-averaged 2008 Suzaku spectrum. 
In all subsequent fits, we use data taken in the 0.6-10.0 and 15.0- 
50.0 keV energy bands for the XIS and HXD respectively. The 
cross-normalisation between the HXD PIN and the XIS detectors 
was accounted for by the addition of a fixed constant component at 
a value of 1 . 1 8 for the HXD nominal pointing position (2005 obser- 
vation) and 1.16 for the XIS nominal pointing position (2008 ob- 
servations); values derived using Suzaku observations of the Crab 
(Ishida, Suzuki & Someya 2007 4 ). We also ignored all data from 
1.7-2.0 keV so as to avoid any contamination from the absorption 
edge due to silicon in the detectors. Since the XIS data were binned 
up at the HWHM of the resolution of the detector, we included a 2 
per cent systematic error in all fits to account for the high statistical 
weight of the bins at low energies. 



6.1 The Broad-Band Suzaku Model 

We firstly applied the best-fitting HETG continuum model (from 
Section 5) consisting of a power-law continuum, a partial-coverer, 
the fully-covering warm absorber (zones 1-4), Fe K emission lines 
and a black body to parameterise the soft excess to the time- 
averaged 2008 broad-band Suzaku data from 0.6-50.0 keV when 
the source was observed to be in a period of high flux (Fo.5-10.0 = 
5.03 x 10" 11 erg cm" 2 s" 1 ). A historical lightcurve of NGC 4051 
in the X-ray band is shown in Figure 9 where it can be seen that the 
2008 observations were made during a period of relatively high flux 
for this source. Since the 2008 data were obtained in two separate 
observations due to scheduling constraints, we initially only con- 
sidered the much longer ~275 ks exposure. We included all of the 
soft X-ray, fully covering absorption zones required by the HETG 
data (including the soft emission lines) and fixed the parameters 
at the best-fitting values listed in Table 3 since the 2008 Chandra 
and Suzaku data were quasi-simultaneous. This resulted in a very 



4 ftp://legacy.gsfc. nasa.gov/suzaku/doc/xrt/suzakumemo-2007-1 l.pdf 
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Figure 9. A historical lightcurve showing the 2-10 keV flux from all ob- 
servations of NGC 4051 with ASCA, BeppoSAX (simultaneously with 
ROSAT), Chandra, XMM-Newton and Suzaku. The lightcurve shows that 
the flux of this source varies significantly over time and that the 2008 Chan- 
dra and Suzaku observations that we describe here (top-right) were made 
during a period of high flux. Note that the 2005 Suzaku data modelled in 
Section 6.5 as well as previously by Terashima et al. (2009) were taken 
during a period of much lower flux. Data for the lightcurve were obtained 
from Guainazzi et al. (1996), Guainazzi et al. (1998), Uttley et al. (1999), 
Collinge et al. (2001), Uttley et al. (2003), Ogle et al. (2004), Steenbragge 
et al. (2009), Terashima et al. (2009), Shu, Yaqoob & Wang (2010) and 
Turner et al. (2010). 

poor fit with x 2 /d.o.f. — 1351/188. Allowing the model to renor- 
malise slightly improved the fit but still returned a poor fit statis- 
tic of x 2 /d-o.f. = 1165/188 with the main contribution to the 
residuals arising from the significant hard excess seen at energies 
>10keV. This hard excess can be seen in Figure 13 (see Section 
6.5) where broad-band Suzaku spectrum is shown (in red) unfolded 
against a power-law with V = 2. 

6.2 The Suzaku Fe K Profile 

In order to accurately analyse the Suzaku Fe K profile, it was there- 
fore important to firstly parameterise the broad-band continuum 
emission using both the XIS and HXD data. To do this, we mod- 
elled the hard excess by including the PEXRAV model (Magdziarz 
& Zdziarski 1995); an additive component incorporating the re- 
flected continuum from a neutral slab. We tied the photon index and 
the normalisation of the unabsorbed power-law continuum to that 
of the power-law component incident upon the reflector and fixed 
the elemental abundances to Solar (Anders & Grevesse 1989). We 
fixed the cosine of the inclination angle of the source to 0.87 (cor- 
responding to 30°) and tied the folding energy to the cutoff energy 
of the power-law at 300 keV, consistent with no cut-off in the HXD 
PIN band. We found that the combination of the PEXRAV and black- 
body models (to account for the hard and soft excess respectively) 
with the partial-coverer from the HETG model (to account for the 
additional spectral curvature) was able to smoothly paramaterise 
the broad-band continuum well. 

We then began to parameterise the Fe K profile by modelling 
the most prominent emission line corresponding to Ka emission 
from near-neutral iron with a Gaussian. This emission line is found 
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Table 4. The best-fitting Fe K parameters from the Chandra HETG. All values are given in the rest frame of the source. See Section 5.2 for details. 
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Figure 10. Plot showing the ratio of the data-to-model residuals of the 
2008 Suzaku XIS data in the Fe K band in the observed frame. The dotted 
lines show the best-fitting centroid energies of the five significant emission 
/ absorption lines required by the data: from left to right, neutral Fe Ka 
emission at 6.41 keV (red), FeXXV forbidden emission at 6.62 keV (blue), 
blueshifted absorption most likely from FeXXV at 6.81 keV (cyan) and 
emission and blueshifted absorption from FeXXVI Lya at 6.97 keV and 
7.12keV respectively (both magenta). 

to have a centroid energy of En ne = 6.41 ± 0.01 keV in the Suzaku 
data and an intrinsic width of a < 53 eV (FWHM <5 700 km s" 1 ). 
It also has an associated equivalent width of EW = 75tg 4 eV 
against the observed continuum. The inclusion of this component 
improves the fit statistic by A^ 2 = 228. The fit was further im- 
proved (Ax 2 = 12) by including a second line to model the for- 
bidden transition from Fe XXV found in Section 5.2 which appears 
to be manifested in the slight blue wing of the Fe Ka line shown 
in Figure 10. The line width was fixed at a = OeV (intrinsically 
narrow) and although the centroid energy was difficult to constrain, 
it was found to have a best-fitting value of Eu nc = 6.621q'q2 keV. 
The best-fitting values of the line flux and associated equivalent 
width (Fii nc = 5.51±2.39 x 10" 6 photons cm -2 s _1 and EW = 
25^1! eV respectively) are also largely consistent with the corre- 
sponding values found in the Chandra HETG data. 

We also included a third line to model the emission from the 
Fe XXVI \s-2p Lya transition. Including an additional line fixed at 
the best-fitting value found in the HETG data of E rinc = 6.97 keV 
with an intrinsically narrow width of a = OeV (see Section 5.2 
and Table 4) is marginally required by the data (Ax ~ 5) and 
corresponds to an equivalent width of EW = 9.9l 7 'g eV against 
the observed continuum. This line appears to be less significant 
in the Suzaku XIS data than in the Chandra HETG data possi- 
bly due to the strength of two significant absorption lines with 



blueshifted centroid energies of ~6.8 keV and ~7. 1 keV (in the ob- 
served frame) possibly masking some of the emission from the ls- 
2p transition of H-like Fe in the XIS data (see Section 6.3). We also 
note that a slight positive residual can be observed at ~7.5 keV (see 
Figure 10) which is likely due to Ka fluorescence from near- neutral 
Ni which has rest-frame centroid energy of ~7.47 keV. However, 
modelling this feature with a Gaussian is statistically unrequired 
by the data. 

Finally, although no neutral Fe Kj3 emission was statistically 
required by the data (A^ 2 ~ 2), we modelled it for consistency 
with the addition of a further Gaussian. We fixed the line energy at 
7.06 keV, tied the intrinsic width to that of the corresponding Fe Ka 
line and fixed the line flux at 13 per cent of the Ka flux, consis- 
tent with the theoretical flux ratio for near-neutral iron (Kaastra & 
Mewe 1993). This line is unrequired by the data but for consistency 
we keep the line modelled in all subsequent fits. 

6.3 Highly Ionised Absorption 

Upon modelling the Fe K emission it can be observed that signifi- 
cant negative residuals are apparent in the data at energies ~7 keV 
suggestive of the presence of highly ionised absorption lines. At a 
first glance, the inferred blueshift of these lines could suggest that 
these lines are the absorption signature of a high velocity outflow, 
such as that detected by Pounds et al. (2004a). A plot of the ratio of 
the data-to-model residuals in the 2008 Suzaku XIS data is shown 
in Figure 10. We initially parameterised these absorption lines with 
Gaussians with negative fluxes. The first line is found to have a 
centroid energy of Bi inc = 6.81+o;q5 keV, improving the fit by 
Ax 2 = 26 (also see Figure 11). Its intrinsic width was fixed at a 
value of (7 = eV. This absorption line has an associated equivalent 
width of EW = —(2St\\) eV and is most likely due to absorption 
from the Fe XXV resonance line at 6.70 keV. If so, the blueshift of 
the line would suggest that the absorbing material is outflowing 
with a velocity of v out = — (5 500)^2 ioo km s" 1 relative to the 
host galaxy, consistent with the value found by Terashima et al. 
(2009) and Pounds & Vaughan (2010). The second line appears to 
correspond to a centroid energy of Ey me = 7.12 ± 0.04 keV with a 
corresponding equivalent width of EW = —(43^^) eV. The ad- 
dition of this line greatly improves the fit statistic by A^ 2 = 56 
(also see Figure 11). Again the intrinsic width of the line was fixed 
at o — eV. If this line is associated with the \s—2p Lya doublet 
from Fe XXVI at ~6.97 keV then it would correspond to an outflow 
velocity of v out = — (75001J 70 q) kms -1 again consistent with 
the findings of Pounds et al. (2004a) and Terashima et al. (2009). 

After parameterising the absorption lines with simple Gaus- 
sians, we attempted to self-consistently model the lines with the 
addition of a photo-ionised grid of absorbing gas using the XS- 
TAR code of Kallman & Bautista (2001). We adopted a turbulence 
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Figure 11. Plot (in the observed frame) showing the significant contribu- 
tion to the x 2 value of the two highly ionised absorption lines at ~6.8 and 
~7.1keV most likely originating from FeXXV and FeXXVI respectively. 
The Fe K emission lines are already modelled, as detailed in Section 6.2. 

velocity of a — 3 000 km s -1 so as not to saturate the absorp- 
tion lines and to explain their high observed equivalent widths. 
We find that the addition of this zone of gas improves the fit 
by Ax 2 = 113 and requires a best-fitting column density of 
Nh = 8.4_ 2 x 10 cm and a best-fitting value for the ionisa- 
tion parameter of log £ = 4. l^o' 2 • The blueshift of the zone corre- 
sponds to an outflow velocity of v out = — (5 800l* 6 2 y ) kms -1 
(~— 0.02c), consistent with the values found by Pounds et al. 
(2004a) and Terashima et al. (2009) for the highly ionised outflow. 
We note that no further zones of absorption are statistically required 
by the data to model the absorption at FeK and that we do not 
see any evidence of the higher velocity (v out ~ — 30 000 km s -1 ) 
component reported by Pounds & Vaughan (2010) (see Figure 12). 
Indeed, by fixing two additional Gaussians centred on energies of 
7.5 and 7.7 keV (corresponding to the observed energies reported 
by Pounds & Vaughan 2010), we note that these are unrequired 
by the data and have upper limits on their equivalent widths of 
EW < 7 eV and EW < 13 eV respectively. 

Upon modelling the Fe K emission and the highly ionised ab- 
sorption, the reflection scaling factor of the PEXRAV component 
was found to have a best-fitting value of R — 1.231q'qJ?. We note 
that a sub-Solar Fe abundance appears to be preferred by the data 
with the value dropping to Af c = 0.37to'5i times Solar. How- 
ever, this could simply be due to the limitations of fitting the spec- 
trum with a simple neutral slab. We also note that the ratio of the 
normalisations of the absorbed to total power-law components at 
1 keV (~0.3) appears to be consistent with the best-fitting values 
obtained from the HETG spectrum (i.e. still suggestive of a ~30 per 
cent covering fraction for the partial-covering zone; also see Sec- 
tion 5.1). The final best-fitting parameters of the broad-band model 
to the Suzaku XIS/HXD data are shown in Table 5. The overall fit 
statistic corresponds to \ 2 /d.o. f. — 218/187. 

Finally, as a consistency check, we returned to the HETG 
data to test for the presence of the two highly ionised absorption 
lines which we detect with Suzaku. Taking the best-fitting HETG 
model from Section 5, we added two Gaussian profiles with neg- 
ative fluxes fixing the centroid energies at En nc — 6.81 keV and 
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Figure 12. Plot showing the FeK region for the three Suzaku observa- 
tions of NGC 405 1 . We find strong evidence for a highly ionised outflow 
with D ou t ~ — 6 000 km s -1 due to the two absorption lines at ~6.8 and 
~7.1 keV which we associate with the ls-2p transitions from Fe XXV and 
FeXXVI respectively. However, we see no evidence for the much faster 
component with Vout ~ — 30 000 km s -1 reported by Pounds & Vaughan 
(2010) which they associate with the presence of two additional absorption 
lines at ~7.5-7.7 keV. We note that slight residuals at ~8.0 and ~8.3 keV 
can be seen in the 2005 data; these are likely due to the ls-3p transitions 
of Fe XXV and Fe XXVI respectively but are not statistically significant. The 
2005 data are shown in black and the two separate 2008 observations are 
shown in red and green corresponding to the 275 and 78 ks exposures re- 
spectively. 

E\i nc = 7.12 keV respectively and the intrinsic width of the two 
lines at a — OeV, consistent with the Suzaku data. Neither of 
the two lines were statistically required by the HETG data but we 
find upper limits on the magnitudes of the equivalent widths of 
EW < 21 eV and EW < 41 eV for the two respective lines, 
indicating that they are consistent with the values that we find with 
Suzaku. 

6.4 Absorbed Reflection Model 

Since the PEXRAV and blackbody components served as simple pa- 
rameterisations of the Compton-scattered reflection hump off neu- 
tral material at energies >10keV and the soft excess at energies 
<2keV respectively, we attempted to model these features using a 
more physical approach. We began by removing the PEXRAV com- 
ponent, the blackbody and the three Gaussians at 6.41, 6.62 and 
6.97 keV (modelling emission from near-neutral FeKa, FeXXV 
and FeXXVI respectively). We then replaced these components 
with two absorbed reflection models (consistent with the full broad- 
band model described in Turner et al. (2010)) using the REFLIONX 
code of Ross & Fabian (2005); one low £ (near-neutral) component 
to model the near-neutral Fe Ka emission line and the hard excess 
and one highly ionised component to account for the spectral curva- 
ture in the soft band (i.e. the soft excess) plus the highly ionised Fe 
emission lines. The REFLIONX code models the emergent spectrum 
from a photo-ionised optically-thick slab of gas when irradiated by 
a power-law spectrum and consists of both the reflected continuum 
and line emission for the astrophysically abundant elements. It as- 
sumes a high-energy exponential cut-off of E cut = 300 keV and 
uses the abundances of Anders & Ebihara (1982). All absorption 
was again modelled using the XSTAR 2. llnl 1 code of Kallman & 
Bautista (2001) with the soft X-ray absorption fixed at the best- 
fitting values from the HETG data (Table 3) and the zone of highly 
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Table 5. Table showing the best-fitting rest-frame parameters of the broad- 
band Suzaku XIS+HXD model described in Section 6. "Primary power- 
law continuum: T, photon index; Normalisation. 6 Absorbed power-law. 
C PEXRAV neutral reflector: R, reflection scaling factor; Ape, iron abun- 
dance with respect to Solar. d Blackbody: fcT, temperature; Normalisa- 
tion defined as L39 /-D 2 where L39 is the source luminosity in units of 
10 39 ergs _1 and £>io is the distance to the source in units of lOkpc. 
e Emission/Absorption lines: i?ii ne , centroid energy; a, intrinsic width; 
F line , line flux; EW, equivalent width, -f Highly ionised absorption: iV H , 
column density; £, ionisation parameter; twt, velocity shift. The symbol I 
signifies that the parameter is linked to all parameters with the correspond- 
ing number. Note that the high £ XSTAR zone replaced the highly ionised 



absorption lines and they were not modelled simultaneously. 
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ionised absorption detected with Suzaku consistent with the values 
obtained in Section 6.3. Any soft X-ray emission lines unable to be 
accounted for by the ionised reflector were modelled with Gaus- 
sians as required. This model can be expressed as: 

F(E) = WABS X (XSTARhetG X XSTARp c K X [PL int . 

+ (PL X XSTARp C ) + (REFLIONXi ow £ X xstar) 
+ (REFLIONX h igh£ X XSTAR)] + GA^'J, (2) 



where WABS is the absorption due to Galactic hydrogen in our 
line-of-sight, XSTARhetg corresponds to the four fully-covering 
soft X-ray absorption zones fixed at the best-fitting values from the 
HETG data (see Table 3), XSTARp c k is the zone of highly ionised 
absorption (see Section 6.3), PLi nt . is the intrinsic ("unabsorbed") 
power-law continuum, PL x XSTAR pc corresponds to the partial- 
coverer, the REFLIONX components are the near-neutral and highly 
ionised absorbed reflectors (each of which is absorbed by a layer of 
gas within XSTAR, without which the soft excess cannot be mod- 
elled; see below) and GA'„' corresponds to the soft emission lines 
which the reflector is unable to account for (also see Section 4.3); 
primarily the forbidden transitions from O VII, Ne IX and Si XIII. 
These lines were fixed at the best-fitting values from the HETG 
data (listed in Table Al). This forms a highly complex model since 
were are attempting to account for a complex ionised scattering 
spectrum with components that cannot neccessarily account for all 
of the physical effects that we expect to be present. Therefore, its 
apparent complexity may simply be an artefact of the limitations of 
the current physical models available, e.g. compared to physically 
realistic disc wind models (see Sim et al. 2008, 2010). 

We fixed the redshift values of the two reflectors to that of 
the host galaxy (i.e. z = 0.002336), tied the iron abundances to- 
gether and tied the photon indexes to that of the incident power- 
law continuum. No additional velocity broadening was applied to 
the reflected spectrum. We find that the photon index steepens to 
T = 2.50 ± 0.02. The ionisation parameter of the low £ reflec- 
tor pegs at a value of £ = 10 erg cm s~ 1 (the lowest value allowed 
by the model) whereas the highly ionised reflector modelling the 
soft excess requires a best-fitting value of log£ = 3. 81^1 . We 
note that the highly ionised reflector is able to account for the 
soft excess due to its enhanced reflectivity in the soft band at high 
levels of ionisation although we stress that the soft excess cannot 
be modelled by the REFLIONX component alone without introduc- 
ing a zone of absorption (in addition to the four warm absorbers, 
the partial-coverer and the Galactic absorption). The absorber in 
front of the reflector is responsible for producing deep bound-free 
edges which reduce the observed reflected flux at higher energies, 
thus revealing positive spectral curvature in the soft band. Such 
a component may be representative of scattering off a disc wind. 
The highly ionised reflector also removes the residuals at Fe K and 
well models the ionised emission lines from Fe XXV and Fe XXVI. 
This model appeared to leave very few residuals in the data and 
gave a good fit with x 2 /d.o.f. = 181/187, a significant improve- 
ment over the phenomenological model described in Section 6.3 
(X 2 /d.o.f. = 218/187). The best-fitting values of the reflection 
components and their associated zones of absorption are shown in 
Table 6. 



6.5 Long-Term Spectral Variability 

In order to assess the nature of the long-term X-ray spectral vari- 
ability of NGC 405 1 , we attempted to simultaneously model the 
2008 data with the 2005 Suzaku data over the broad-band 0.6- 
50.0 keV energy range. The 2005 data have previously been de- 
scribed in detail by Terashima et al. (2009) where NGC 4051 
was found to be in the low-flux state with a 0.5-10.0 keV flux of 
F = 1.32 x 10~ n ergem -2 s _1 compared to the much higher 



flux of F = 5.03 x 10- n ergcm- 2 s 



1 in 2008. The spectrum 
was observed to be much harder in the low-flux state with a signifi- 
cant excess at energies > 10 keV. Figure 13 shows the relative fluxes 
of the three observations. The large changes in spectral shape can 
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Table 6. Table showing the best-fitting rest-frame parameters of the broad- 
band Suzaku XIS+HXD absorbed reflection model described in Sec- 
tion 6.4. "Primary power-law continuum: T, photon index; Normalisa- 
tion. b Absorbed power-law: ./Vh, column density; ?, ionisation parameter. 
c Neutral reflector: Ap c , iron abundance with respect to Solar. d Zone of ab- 
sorption. The redshift of the zone was fixed at the same value as that of the 
host galaxy (i.e. z = 0.002336). Tonised reflector. The symbol p signifies 
that the parameter has pegged at the maximum / minimum value allowed 
by the model and the symbol I signifies that the parameter is linked to all 
parameters with the corresponding number. 
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be well described by a two-component model consisting of a soft 
variable component superimposed over a hard constant component 
(Miller et al. 2010). 

We simultaneously applied the best-fitting model described 
in Section 6.4 to the 2005 and 2008 XIS and HXD Suzaku data. 
We also included the second, shorter exposure (78 ks) observa- 
tion from 2008 in our joint-analysis. We tested to see if the fully- 
covering warm absorber parameters varied between observations 
but the best-fitting values were consistent within the errors across 
the 2005 and 2008 data. We tied together all model parameters 
between the three observations except for the relative normalisa- 
tions of the absorbed and intrinsic power-law components and the 
near-neutral and ionised reflection components (REFLIONX) which 
we allowed to vary to account for the long-term spectral variabil- 
ity. We find that the normalisations of the absorbed power-law and 
the near-neutral reflection component appear to be largely constant 
across the three observations (see Table A2), consistent with the 
findings of Miller et al. (2010). We note that slight changes in the 
flux of the highly ionised reflector can be observed between 2005 
and 2008. However, we find that the spectral variability can ulti- 
mately be described by large changes in the normalisation of the 
intrinsic unabsorbed power-law component which increases by a 
factor of ~7 from the low-flux state in 2005 to the high-flux state 
in 2008. This corresponds to a high covering fraction of the partial- 
coverer of ~70 per cent in the 2005 low-flux data compared to a 
much lower covering fraction of ~30 per cent when in the high- 
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Figure 13. A plot showing the relative fluxes of the combined Suzaku XIS 
and HXD PIN data for NGC 4051 in the observed frame. The data are un- 
folded against a power-law with T = 2. The 2005 data are shown in black. 
The two separate 2008 observations are shown in red and green correspond- 
ing to the 275 and 78 ks exposures respectively. 

flux state and supports the notion of an empirical two-component 
model whereby the constant reflection component contributes to the 
majority of the hard X-ray flux whereas changes in the normalisa- 
tion of the intrinsic power-law account for the variations in spectral 
shape at lower energies. The 0.5-100 keV fluxes corresponding to 
the relative normalisations of the power-law and reflection compo- 
nents are given in Table 7. A plot of the relative contributions of the 
individual model components to the joint-fit is shown in Figure 14. 

We note that upon considering the joint-fit of all Suzaku ob- 
servations, we also find that an additional emission line at i5n nc = 
5.44 ± 0.03 keV is statistically required in the 2005 data, possibly 
attributed to CrKa. Modelling this with a Guassian profile with an 
intrinsically narrow fixed width of a = e V provides a measure- 
ment of the equivalent width of EW = 42l 12 eV (consistent with 
Turner et al. 2010) and improves the fit statistic by A\ 2 = 20. We 
also include this line in the 2008 data and find that the flux of the 
line appears to be consistent with the 2005 data with the equivalent 
width dropping to EW = 15_g eV against the observed contin- 
uum as expected. This suggests that the emission line may have 
remained constant in flux over this three-year period. A detailed 
analysis of the nature of this line is given by Turner et al. (2010) 
with its possible origins discussed in a companion paper (Turner & 
Miller 2010). 

This is then considered to be our final model which appears 
to describe the data well with all remaining free parameters across 
all three observations consistent with the values listed in Table 6 
resulting in a final fit statistic of \ 2 /d.o.f. = 560/549. 



7 DISCUSSION 

7.1 Highly Ionised Outflow 

The Suzaku XIS spectrum of NGC 4051 reveals two absorption 
lines with centroid energies at ~6.8 and ~7.1keV, as described 
in Section 6.3. This suggests that the absorption is significantly 
blueshifted and could be the signature of a highly ionised, high 
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Figure 14. A plot showing the relative contributions of the individual model 
components to the three broad-band 0.6-50.0 keV Suzaku spectra in the 
observed frame. The data are shown in black with the sum of all model 
components superimposed in red. The absorbed and unabsorbed power-law 
components are shown in cyan and navy blue respectively. The absorbed 
neutral and ionised reflectors are shown in magenta and green respectively. 
Note that the absorbed power-law component and the reflection components 
appeal' to remain largely constant in flux across the three observations with 
the changes in spectral shape instead accounted for by significant variations 
in the normalisation of the intrinsic unabsorbed power-law component. 



Table 7. The upper-half of the table shows the fluxes corresponding to 
the best-fitting values for the power-law (both absorbed and unabsorbed), 
reflected continua and FeKa components obtained in Section 6.5. The 
fluxes are given across the 0.5-100 keV energy band and have units 
10~ 12 erg cm -2 s -1 . The covering fraction of the partial-coverer is also 
shown for each observation. The lower-half of the table shows these values 
expressed as a ratio compared to the values obtained for the 2005 Suzaku 
data. Obs. 1 corresponds to the 2005 data whereas Obs. 2 and 3 correspond 
to the two 2008 observations (275 and 78 ks respectively). It can be seen 
that the fluxes of the reflection components, the absorbed power-law and 
the neutral FeKa emission line are consistent with remaining largely con- 
stant across all epochs whereas large variations in the flux of the intrinsic 
unabsorbed power-law are apparent between observations with the flux con- 
siderably lower in 2005 when the source was in a historically low-flux state. 
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velocity outflow. Modelling with an XSTAR grid with a turbulence 
velocity of a = 3 000 km s" 1 requires one zone of absorption with 
a best-fitting column density of Nn = 8. 4_ 2 x 10 cm and an 
ionisation parameter of log£ = 4.1+°'^ (as per Section 6.3). As- 
suming that the lines correspond to K-shell absorption of Fe XXV 
and FeXXVI, the blueshift of the zone corresponds to an outflow 
velocity of v out — — (5 800+^200) kms -1 , largely consistent with 
the findings of Pounds et al. (2004a) and Terashima et al. (2009). 
It is also conceivable that this is the highly ionised signature of the 
high velocity outflow detected by Steenbrugge et al. (2009) with 
the Chandra LETG which they deduce to have a comparable out- 
flow velocity of v out ~ — 4 500 km s _1 although with an ionisation 
parameter ~10 times lower (i.e. log £ ~ 3.1). 



7. 1. 1 The Kinematics of the Highly Ionised Absorption 

For a homogeneous, radial, outflowing wind, the column density 
along the line of sight is given by A^h = J^° ut n e (R)dR, where 
Rin and Rout are the inner and outer radii along the line of sight 
through the wind respectively. Assuming we are looking down a 
spherically symmetric homogeneous wind towards the inner radius, 
then Rout tends to 00. So, combining this with the definition of 
the ionisation parameter (equation 1) leads to Ri n = Li on /Nn£,- 
However, if the thickness of the wind is defined as AR = R ou t — 
Ri n then the inner radius may extend inwards if the wind is clumpy 
(i.e. if R » AR). This yields: 



Rin < 



Lie 



(3) 



A value for the ionising luminosity can be calculated by 
extrapolating the total power-law continuum of the broad-band 
Suzaku model from the 2008 data and integrating from 1 to 1 000 
Rydbergs. This gives a value of L lon ~ 1.4 x 10 43 ergs -1 . Com- 
bining this with the best-fitting values of the column density and the 
ionisation parameter from the XSTAR model (Table 5) then yields 
a value for the radius of the absorbing material of R^ n < 1.3 x 
10 16 cm which is of the order of < 0.004 pc or < 5 1-d from the 
central engine. Assuming that the material escapes to infinity then 
a lower limit on the radius can also be calculated by considering the 

outflowing material's escape velocity where v eBC = 2G J >I . Rear- 
ranging for R gives a lower limit of R > 1.4 x 10 15 cm placing 
the constraints on the radius of the highly ionised absorber to be 
0.5 < R < 5 1-d from the central source and infers a constraint on 
the electron density of 7 x 10 6 < n c < 6 x 10 s cm -3 (equation 
1). The constraints on the radius derived here appear to be consis- 
tent with the radius of 0.5-1 1-d found by Krongold et al. (2007) for 
their highly ionised zone of absorption found in the XMM-Newton 
RGS spectrum. This is well within the 10-15 1-d dust sublimation 
radius in NGC 4051 (see Krongold et al. 2004 fig. 10) and so can- 
not have an origin in the molecular torus. Instead, as the location of 
the highly ionised absorber also appears to be contained within the 
location of the He II broad emission line region (BELR; Krongold 
et al. 2007), an accretion disc wind is perhaps a more likely origin 
for this absorption signature. 

By considering the case of a homogeneous, spherical flow un- 
der the assumption that the outflow velocity remains approximately 
constant on the compact scales observed here (although decelerat- 
ing at some larger radius), an estimate on the mass outflow rate of 
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the highly ionised zone can be calculated from the simple conser- 
vation of mass: 

Mout = 47r6^Y^Woutm p , (4) 

where Lion/£ = n c R 2 . Here, M ou t is the mass outflow rate 
in units g s~ subtending a solid angle, Airb (where < b < 1), 
where L- lon is is the ionising luminosity from 1 to 1 000 Rydbergs in 
units ergs -1 , £ is the ionisation parameter in units erg cms -1 , v out 
is the outflow velocity in units cms -1 and m p is the rest mass of 
a proton. Initially assuming that the absorber is fully covering (i.e. 
b = 1) and taking the best-fitting values of the ionisation parameter 
and outflow velocity given in Table 5 and the value for the ionising 
luminosity derived above yields an estimate on the mass outflow 
rate of M out ~ O.2M yr -1 . This value then translates into an 
estimate of the kinetic output of the outflow via: 

1 • 2 

-Bout = -M out v out , (5) 

which corresponds to a value of E out ~ 2 x 10 42 ergs -1 , 
which is approximately 3-4 per cent of the bolometric output of 
the AGN (Ogle et al. 2004). The momentum rate of the outflow is 
then Moutvout ~ 8 x 10 33 gems -2 . 

For a momentum-driven outflow whereby r ~ 1 (i.e. where 
the photons scatter once before escaping and so all momentum 
is transferred to the outflow) accreting at near-Eddington (i.e. 
Mout ~ Medd), the total wind momentum flux must be of the 
same order as the photon momentum flux (King 2010), i.e.: 

»> Bedd 

M ut«out ~ • (6) 

c 

Assuming a black hole mass of M B h = 1-73 x 10 6 M 
(Denney et al. 2009) for NGC 4051, the Eddington luminosity is 
found to be Ledd ~ 2.2 x 10 44 ergs -1 which yields Ledd/c ~ 
7.5 x 10 33 erg cm -1 ; a value comparable to the calculated outflow 
momentum rate of MoutWout ~ 8 x 10 33 gems -2 . So it can be 
seen that equation 6 is largely satisfied if we assume that the highly 
ionised outflow is momentum-driven by Thomson scattering such 
as in the similar case of the quasi-stellar object (QSO) PG 121 1+143 
(King & Pounds 2003; Pounds et al. 2003). 

Furthermore, the Eddington luminosity is also defined by: 

£edd = »?Meddc 2 , (7) 

where r\ is the efficiency of mass-to-energy conversion and 
takes a value between and 1. Combining this with equation 6 
yields an expression for the wind velocity (King et al. 2010): 

V 

Vout ~ — c, (8) 

m 

where rh = Mjcc/Medd is the accretion rate in Eddington 
units. Since v out = 0.02c here, a value of r\ ~ 0.02 would be re- 
quired to sustain the Eddington luminosity of this source if it were 
accreting at the Eddington rate (i.e. rh ~ 1). Note that the max- 
imum value of 7] is 0.06 for a non-rotating, Schwarzschild black 
hole. 



The bolometric luminosity of NGC 405 1 calculated by Ogle et 
al. (2004) and corrected for the distance of 15.2 Mpc assumed here 
is Lboi ~ 7 x 10 43 ergs -1 . Comparing this to the Eddington lu- 
minosity suggests that this object has an accretion rate of only ~30 
per cent of Medd which would correspond to an observed mass 
accretion rate of M acc = L bol /r/c 2 ~ 0.06 M Q yr -1 assuming an 
efficiency of r\ — 0.02. However, the observed luminosity of the 
system with respect to Eddington may be underestimated if some 
of the emission is obscured by a Compton-thick wind and scattered 
out of the line of sight (a notion perhaps supported by the appar- 
ent presence of high column densities towards this source). So it is 
perhaps conceivable that NGC 405 1 is accreting at an appreciable 
fraction of the Eddington rate. Furthermore, Sim et al. (2010) show 
that radiatively-driven accretion disc winds can exist with similar 
properties to those that we observe here for sources with luminosi- 
ties of a few tens of per cent of Ledd- 

7.1.2 The Covering Fraction 

From the consideration of mass conservation (equation 4), it can be 
seen that the mass outflow rate is dependent upon the covering fac- 
tor, b, which is likely to be some significant fraction of 4-7T sr. The 
calculations in Section 7.1.1 were based on the assumption that the 
outflow is a radial flow in a spherical geometry with 6=1. How- 
ever, if the outflow takes the form of a bi-conical geometry (Elvis 
2000), as proposed by Krongold et al. (2007) for this source then 
the covering fraction may be lower. 

The Chandra HETG and Suzaku XIS spectra reveal the pres- 
ence of two highly ionised Fe K emission lines most likely originat- 
ing from Fe XXV and Fe XXVI. If these lines are the signature of the 
same ions associated with the highly ionised, high velocity outflow, 
we can attempt to obtain an estimate on the covering fraction of the 
material. Bianchi & Matt (2002) calculate the predicted equivalent 
widths of the Fe XXV and Fe XXVI emission lines against the total 
(primary + reprocessed) continuum for a given photon index and 
column density. Figures 8 and 9 of Bianchi & Matt (2002) show 
that for a column density of Afj = 10 23 cm -2 and a photon in- 
dex of r = 2.4, the expected equivalent widths are ~50eV and 
~15eV for the FeXXV and FeXXVI emission lines respectively. 
As the column density of the highly ionised absorption that we ob- 
serve here with Suzaku is Nh ~ 8 x 10 22 cm -2 (see Section 6.3) 
and the best-fitting photon index is T ~ 2.5 (see Section 6.4) in the 
2008 data, we can make a direct comparison. Within the errors, the 
equivalent widths of ~25 eV and ~ 10 eV obtained for the Fe XXV 
and Fe XXVI emission in the Suzaku data respectively (see Table 5) 
suggest that the values are on the order of ~50 per cent lower than 
the predicted maximum values perhaps implying that the covering 
factor is on the order of b ~ 0.5. 

Furthermore, Tombesi et al. (2010) analysed a sample of 42 
radio-quiet AGN observed with XMM-Newton. They detected 22 
absorption lines at rest-frame energies >7. 1 keV, implying that high 
velocity outflows may be a common phenomenon in radio-quiet 
AGN. The fraction of AGN having at least one feature with a blue- 
shifted velocity in their sample is 17/42, corresponding to ~40 per 
cent of the objects. They also find that the global covering fraction 
of the absorbers that they detect is estimated to be in the range C ~ 
0.4-0.6, suggesting that the outflowing winds generally have large 
opening angles. 

We can also obtain a rough estimate on b by comparing the 
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observed accretion rate of the source (M acc ~ 0.06 M Q yr _1 ; see 
Section 7.1.1) with the calculated spherical outflow rate from equa- 
tion 4 (Mout ~ 0.2 Mq yr _1 ) under the assumption that the mass 
outflow rate either does either not exceed or is of the same order as 
the mass accretion rate (i.e. M>ut < Macc). This suggests a cover- 
ing fraction of b < 0.3 although this could be higher if the observed 
luminosity with respect to Eddington is underestimated somewhat 
perhaps due to partially-covering Compton-thick clouds in our line 
of sight. So it seems that although the absorption component is un- 
likely to be fully covering, the observational evidence suggests that 
the outflow is not highly collimated with b still likely to correspond 
to a significant fraction of 1 which appears to be consistent with the 
detailed reverberation analysis of Miller et al. (2010) who state that 
a global covering fraction of > 40 per cent is required to match the 
observations. 



7.1.3 Implications for Feedback 

Given the high accretion rates required for the growth of super- 
massive black holes, significant outflows may be a consequence 
of near-Eddington accretion (King 2010). If so, they could be an 
important factor in galaxy evolution and in establishing the ob- 
served M-o and M-Mbui gc relations between the supermassive 
black hole and its host galaxy (Magorrian et al. 1998; Ferrarese & 
Merritt 2000; Gebhardt 2000). 

Taking the bolometric luminosity of Lboi ~ 7 x 10 43 erg s _1 
of NGC 4051, we can calculate an estimate for the mass accre- 
tion rate of M acc ~ 0.06 M Q yr _1 . Thus, for the black hole to 
accrete a mass of M BH = 1.73 x 10 6 M (Denney et al. 2009), 
the Salpeter e-folding time would be ~3 x 10 7 yr. As calculated 
in Section 7.1.1, the energy output of the high velocity outflow is 
Sout ~ 2 x 10 42 ergs _1 assuming that b ~ 1. Integrating this 
over the e-folding time of the AGN gives a total energy output 
of -Etot ~ 10 57 erg (even if we assume a conservative value of 
b = 0.1, this still provides a total output of Etot ~ 10 56 erg). This 
value can then be compared with the binding energy of the bulge of 
the host galaxy. Marconi & Hunt (2003) and Haring & Rix (2004) 
provide an observable relationship between the mass of the central 
object and the mass of the galaxy bulge, i.e. Mbh ~ 10~ 3 Mb u i ge . 
This can then be used to calculate the binding energy of NGC 405 1 
whereby B.E. ~ cr 2 M bulgc . Nelson & Whittle (1995) provide a 
value for the velocity dispersion of NGC 4051 of a = 88kms~ x 
which then yields an estimate on the binding energy of B.E. ~ 
3 x 10 56 erg, comparable with the likely total energy deposited by 
the wind (assuming that the wind persists for the lifetime of the 
AGN). Therefore, it can be seen that the total energy output of the 
high velocity outflow could have a considerable influence on its 
host galaxy environment. Indeed it may be conceivable that NGC 
4051 is a low mass analague of the QSOs PG121 1 + 143 (Pounds & 
Reeves 2009) and PDS 456 (Reeves et al. 2009) which may deposit 
as much as 10 60 erg throughout the lifetime of the AGN. 

7.2 The Soft X-ray Absorption / Emission 

The HETG spectrum shows evidence for a wealth of absorption 
lines in the soft X-ray band. Modelling this warm absorber with 
XSTAR requires four distinct ionisation zones of gas with the ion- 
isation parameter ranging from log£ = —0.86 to log £ = 2.97 
and column densities on the order of 10 20 -10 21 cm -2 (see Table 
3). The zones appear to be outflowing with velocities on the order 



of a few 100 km s -1 with the general trend appearing to be that 
the larger outflow velocities appear to correspond to more highly 
ionised zones. This is perhaps suggestive of a geometry whereby 
all of the zones form part of the same extended outflowing wind 
with the higher velocity, more highly ionised components originat- 
ing much closer to the central source. Indeed, as the mass outflow 
rate (equation 4) and the ionisation parameter (equation 1) combine 
to give: 

M out oc ^ lon ^° ut — const. (9) 

(King et al. 2010), such a trend between ionisation parame- 
ter and outflow velocity would be expected to be observed if the 
mass outflow rate is to be conserved. Such a correlation has pre- 
viously been noted in NGC 4051 by Pounds et al. (2004a) and, if 
confirmed, could provide strong evidence of a cooling shock (King 
et al. 2010). 

Using equation 3, we are able to obtain upper limits for the 
radius of the absorbing zones of gas (see Table 8). However, these 
values are largely unconstrained and result in very conservative es- 
timates on the order of kpc for the two most lowly ionised zones. 
The constraints on the medium and higher ionisation zones are a 
little tighter with the radius of the material falling within ~70 and 
~2pc from the central source respectively. We note that the ab- 
sorbers could exist at smaller radii still if the material is somewhat 
clumpy or filamentary. 

Calculations of the mass outflow rate of zones 1-4 using equa- 
tion 4 prove to significantly exceed the Eddington accretion rate 
of the source by several orders of magnitude since they assume 
that the absorber forms part of a fully-covering, homogeneous, ra- 
dial flow. However, if we assume that the mass outflow rate does 
not exceed the observed mass accretion rate since the energy re- 
quired to accelerate the outflow must come from the Mailing gas 
(i.e. Mout < 0.06 M Q yr -1 ), we can scale the soft zones accord- 
ingly and place upper limits on the covering fraction. This results 
in values ranging from b < 4.4 x 10~ 4 to b < 0.56 for the lowest 
and highest ionisation zones in the HETG data respectively. The 
values for the lower ionisation zones are generally consistent with 
the opening angles of the warm absorber components detected by 
Steenbrugge et al. (2009) with the Chandra LETG. We also calcu- 
lated upper limits for the kinetic output of the outflow using our 
assumption about the mass outflow rate stated above and equation 
5. This resulted in values ranging from i5 ou t < 1.8 x 10 39 erg s -1 
to -Eout i$ 2.9 x 10 40 ergs -1 for the five zones of absorption, i.e. 
very small fractions of the total bolometric output of the AGN in 
constrast to the high velocity outflow that we observe with Suzaku 
which has a kinetic output at least two orders of magnitude higher 
and could potentially be highly significant in terms of galactic feed- 
back. These values are summarised in Table 8. 

We also detect several narrow emission features at energies 
<2keV. The ionised reflector is unable to model all of the emission 
lines and, in particular, cannot account for the forbidden transitions 
from O VII, Ne IX and Si XIII. Previous studies with the RGS on- 
board XMM-Newton have claimed the presence of RRC (Ogle et 
al. 2004; Pounds et al. 2004a) associated with H-like and He-like 
ions of elements such as O and Ne, especially when the source 
is found to be in an extended period of low flux, suggesting that 
the soft X-ray emission lines may have a photo-ionised origin. In- 
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Table 8. Table showing the upper limits of the radius, R, kinetic energy out- 
put, _Bk, and covering fraction, b, obtained for the five individual zones of 
absorption detected in the Chandra HETG data (see Section 5.1) assuming 
that the mass outflow rates for each zone does not exceed or is on the order 
of ~0. 1 Mq yr~ 1 . The absorption zones directly correspond to those listed 
in Table 4. We also include the highly ionised zone of absorption detected 
with Suzaku for comparison. 



Absorption 


R 




b 


Component 


(PC) 


(ergs" 1 ) 




Zone 1 


< 8.7 x 10 4 


< 1.8 X 10 39 


< 4.4 x 10~ 4 


Zone 2 


< 1.2 x 10 4 


< 7.4 x 10 39 


< 5.2 x 10~ 3 


Zone 3a 


<75 


< 1.7 x 10 40 


$ o.io 


Zone 3b 


<65 


< 3.8 x 10 40 


<0.05 


Zone 4 


< 1-8 


< 2.9 x 10 40 


<0.56 


High 5 


0.0004 < R < 0.004 


~2 x 10 42 


~0.5 



deed an origin in photo-ionised gas (either formed in the NLR or 
in an extended photo-ionised outflow) has been claimed for soft 
X-ray lines in many other objects (e.g. Bianchi et al. 2006). If the 
emission lines detected here do indeed originate in a photo-ionised 
plasma, an estimate on the electron density of the O VII emission 
(the strongest soft X-ray emission line detected here) can be cal- 
culated using the R ratio for He-like ions of Porquet & Dubau 
(2000), which is defined as R(n e ) = where z is the inten- 
sity of the forbidden line and x and y are the intensities of the 
intercombination lines. An upper limit on the flux of the O VII 
intercombination line(s) can be found by including an additional 
Gaussian in the HETG spectrum with the centroid energy fixed 
at 568.5 eV and the intrinsic width of the line tied to that of the 
associated forbidden transition. This component is unrequired by 
the data and provides an upper limit on the line flux of F Un c < 
2.85 x 10~ 5 photons cm~ 2 s _1 . Combining this with the flux of the 
O VII forbidden line of Fn nc = 1.251^32 x 10~ 4 photons cm s -1 
(see Table Al) provides a lower limit on the value of the R ratio 5 of 
R > 3. Then, by considering the relationship between the R ratio 
and electron density for O VII shown in fig.8 of Porquet & Dubau 
(2000), an upper limit on the electron density for the O VII emis- 
sion can be obtained which is on the order of n c < 10 11 cm - ' . If 
the majority of O VII absorption is being modelled by zone 2 (see 
Table 3 and Figure 7), then by assuming that the emission has an 
ionisation parameter of log £ ~ 0.6, equation 1 allows us to place 
a lower limit on the radius of this emission of R > 3.8 x 10 15 cm, 
which is on the order of ~10 4 R g or ~ 1.5 1-d. This radius is largely 
consistent with the value of R < 100 1-d found by Steenbrugge et 
al. (2009) and the value of R < 3.5 1-d for the lowly ionised com- 
ponent detected by Krongold et al. (2007). 



7.3 The Nature of the Long-Term Spectral Variability 

NGC 405 1 is found to be highly variable on both short- (ks) and 
long-term (years) time-scales with the spectral trend being that the 
X-ray spectrum flattens above a few keV as the source flux drops. 
Occasionally NGC 405 1 is found to fall into an extended period of 
low flux (e.g. Uttley et al. 1999) revealing a very hard X-ray spec- 
trum. The general consensus is that the primary power-law disap- 



5 One caveat to this is that densities obtained with the R ratio are also 
dependant upon photoexcitation of the associated resonance line and optical 
depth effects. 



pears from view as the flux drops leaving the remaining emission 
dominated by a hard reflection component (e.g. Guainazzi et al. 
1998; Pounds et al. 2004a). Such a low-flux state was observed in 
2005 by Terashima et al. (2009) when a very strong hard excess 
was observed at energies >10keV which, if modelled with neutral 
reflection (i.e. a PEXRAV component), returns a value of R ~ 7 
(where R = 1 corresponds to reflection from a semi-infinite slab). 
This is also apparent in Table 7 where it can be seen that the 0.5- 
100 keV flux of the reflector is significantly higher than that of the 
intrinsic power-law in the 2005 data. However, the excess could 
not be accounted for by reflection alone due to the relatively weak 
FeKa emission component (EW ~ 140 eV). Instead, Terashima 
et al. (2009) were able to model the broad-band spectra by intro- 
ducing a partially-covered power-law component which gives rise 
to the hard spectral shape. 

In Section 6.4, we successfully modelled the broad-band 2008 
Suzaku data with a model consisting of an intrinsic power-law 
component, a partial-coverer, near-neutral and ionised reflection, 
ionised absorption from a variety of ionisation states and several 
soft X-ray photo-ionised emission lines. This model accounts for 
the spectral curvature well with the partial-coverer appearing to 
have a low covering fraction of ~30 per cent. We then attempted to 
simultaneously fit all three broad-band Suzaku spectra from 2005 
and 2008 using the same model (see Section 6.5). We find that 
the long-term spectral variability can be accounted for largely by 
changes in the normalisation of the intrinsic (unabsorbed) power- 
law. Since the absorbed power-law remains unchanged, the spectral 
variability cannot be accounted for by simple changes in the cover- 
ing fraction of the partial-coverer alone. If this were the case, then 
the total flux (intrinsic + absorbed) would not vary. However, since 
direct changes in the covering fraction are observed between obser- 
vations, the spectral variability cannot be accounted for by intrinsic 
changes in the luminosity of the source alone either. Therefore, ei- 
ther the covering fraction of the Compton-thin partial-coverer and 
the luminosity of the source are inversely correlated or the changes 
in L are not intrinsic. 

If the covering fraction and luminosity are inversely corre- 
lated, one possible cause of this could be that the hot corona above 
the disc may be varying in area. This could have the effect of con- 
centrating the continuum emission from a smaller region, therefore 
increasing the proportion which is absorbed by the partial-coverer 
and lowering the luminosity, such as is observed in the 2005 data. 
Alternatively, the accretion disc itself may become truncated which 
could also have the effect of lowering the intrinsic luminosity and 
concentrating the continuum emission over a smaller area. Taking 
values of M = 0.06 M yr~\ M BH = 1.73 x 10 6 M Q and a - 1, 
from the Shakura & Sunyaev (1973) solution for a thin accretion 
disc, an estimate can be made on the viscous timescale at the ra- 
dius of the innermost stable circular orbit for a Schwarzschild black 
hole (i.e. 6R g ) for NGC 4051 of t viBC ~ 10 7 s (Frank, King & 
Raine 2002). This suggests that the disc could vary on timescales 
as short as a few months or years and so could be a possible cause 
of the long-term variations which we observe. A third possibility 
remains that if the intrinsic ionising luminosity increases, the ab- 
sorbing clouds may become ionised and so more transparent to the 
continuum emission. This could result in an apparent change in the 
covering fraction of the absorber versus luminosity and therefore 
lead to the observed spectral changes between 2005 and 2008. 

Alternatively, if the observed changes in L are not intrinsic, 
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then another possible reason for the observed changes in the flux 
of the intrinsic power law could be due to the intrinsic continuum 
emission disappearing from view as the covering fraction of a fur- 
ther Compton-thick, variable partial-coverer increases therefore re- 
sulting in an apparent drop in continuum flux. Thus, this leaves be- 
hind a constant, hard, reflection component which then dominates 
the spectrum at low fluxes leading to the spectral trend observed 
by Guainazzi et al. (1998) and Pounds et al. (2004a). If the pro- 
posed Compton-thick clouds form part of the same system as the 
more distant Compton-thin clouds, this could explain the correla- 
tion with the line-of-sight covering fraction of the Compton-thin 
partial-coverer (i.e. as / cov of the Compton-thick clouds increases, 
/cov of the Compton-thin clouds increases too). Such an interpreta- 
tion of accounting for spectral variability by allowing for changes 
in the covering fraction of an absorbing layer of gas has also been 
used as a solution to several other well-studied AGN (e.g. 1H 0419- 
577, Pounds et al. 2004b; NGC 1365, Risaliti et al. 2007; MCG- 
6-30-15, Miller, Turner & Reeves 2008; NGC 3516, Turner et al. 
2008). Then, if the origin of the distant, near-neutral reflector that 
we observe is located outside of the partial-coverer (for instance, 
the molecular torus), the apparent value of R > 1 observed in the 
2005 data could simply be due to a delay in the reflection com- 
ponent responding to a previous apparent brighter state since the 
continuum emission from the central engine may also be obscured 
by the Compton-thick clouds at the distance of the torus. 

These results are consistent with the findings of Miller et al. 
(2010) who performed PCA on the three Suzaku observations of 
NGC 4051 from 2005 and 2008. They were able to decompose the 
spectra into the principle modes of variation and confirmed that a 
constant, hard component is present in the spectra at all times with 
a highly variable soft X-ray component superimposed. This was at- 
tributed to variations in the covering fraction of a layer of absorbing 
gas with a covering fraction of > 40 per cent obscuring the central 
source on long time-scales. They suggested that the strong variabil- 
ity on short (ks) time-scales however, is intrinsic to the source. 

The flux of the narrow, unresolved Fe Ka component appears 
to remain constant between the 2005 and 2008 Suzaku observa- 
tions (see Table 7) with the equivalent width of the line decreasing 
against the observed continuum as the flux level of the continuum 
increases (see Table A2). The constancy of the flux over long-term 
time-scales suggests that near-neutral Fe Ka emission does there- 
fore not respond to apparent changes in the continuum level and is 
perhaps supportive of an origin in very distant material. The con- 
stant flux of the neutral FeKa emission line is also observed in 
the PCA of Miller et al. (2010; see figure 2) despite large apparent 
changes in the continuum flux. 

In addition to the near-neutral FeKa emission from distant 
material, we also find evidence for a weak, broader component 
to the Ka emission in the HETG spectrum, as described in Sec- 
tion 5.2. Taking the FWHM of 16 000ll ° ° ° km s" 1 (i.e. a ~ 
8 000 km s _1 ) of this broad Fe Ka component from the HETG data 
allows an estimate on the radius from the central engine of the emit- 
ting material to be calculated assuming simple Keplerian motion. 
This yields an estimate on the radius of R ~ 1 — 2 x 10 3 R g where 
R g = GM/c 2 . Taking the value of 1.73 x 10 6 M Q for the mass of 
the black hole calculated by Denney et al. (2009), this corresponds 
to a value of R ~ 0.1 1-d from the central source. Again, this 
is largely consistent with the findings of Miller et al. (2010) who 
find evidence of a weak variable hard excess at energies >20keV 



with an associated weak, moderately broadened Fe Ka component 
which appears to respond to changes in the continuum on a 20 ks 
time-scale. This may be related to the ionised reflector which is 
seen to vary somewhat in the Suzaku data (see Table 7) and may be 
associated with reflection occuring off the outer disc or off a disc 
wind similar to that calculated by Sim et al. (2010). Miller et al. 
(2010) also find evidence for time-lags between the hard and soft 
bands which are well described by reverberation from material a 
few light-hours away from the illuminating source. This confirms 
the idea that NGC 405 1 has a substantial amount of optically-thick 
material within a few 100 R s from the central source and is again 
supportive of a partial-covering solution to account for the long- 
term spectral variability. 

Finally, we also note that whether parameterising the broad- 
band Suzaku data with a PEXRAV or REFLIONX model, the data 
seem to prefer a sub-Solar abundance of iron with a best-fitting 
value of Afc ~ 0.3 times Solar. Aside from an intrinsically sub- 
Solar abundance of iron in NGC 4051, one other possible in- 
terpretation could be that this is perhaps indicative of spallation 
(Skibo 1997); a phenomenon whereby Fe nuclei are fragmented 
into lighter nuclei by the impact of high energy cosmic rays, thus 
enhancing the abundances of lower Z atoms such as Cr and Mn. In- 
deed, in Section 6.5 we note that we detect an additional emission 
line at -Bii nc = 5.44 keV, whose centroid energy appears to coin- 
cide with the expected line energy of the neutral Cr Ka transition 
thus suggesting an enhanced abundance of Cr. The presence of this 
line is discussed in a companion paper (Turner et al. 2010) where 
they show that the observed line at ~5.44keV cannot be caused 
by a simple statistical fluctuation. Furthermore, the possible origin 
of this line is discussed further by Turner & Miller (2010) where 
they consider an origin in a transient hotspot but ultimately favour 
a spallation interpretation. 



8 CONCLUSIONS 

Through a detailed analysis of the X-ray spectrum of NGC 405 1 
with the Chandra HETG and the SuzaJcu XIS and HXD instru- 
ments, we are able to fully parameterise the warm absorber find- 
ing zones of absorption ranging in ionisation parameter from 
l°g£ — —0.86 to log ^ = 4.1 and ranging in column density from 
iV H ~ 10 20 cm~ 2 to 7Y H ~ 10 23 cirT 2 . The soft X-ray absorber 
zones appear to be outflowing with velocities on the order of a few 
100 km s -1 with the general trend being that the outflow velocity 
increases with ionisation parameter. The kinetic output of the zones 
then appear to be on the order of -E ou t < 10 39 -10 40 ergs -1 ; val- 
ues which correspond to very small fractions of the total bolometric 
output of the AGN and are negligible with respect to the binding en- 
ergy of the galactic bulge (although such kinetic luminosities could 
still significantly influence the ISM; Hopkins & Elvis 2010). 

Regarding the 2008 Suzaku spectrum of NGC 4051, we de- 
tect the presence of two statistically significant absorption lines at 
~6.8 and ~7.1 keV. If these lines are the signature of He-like and 
H-like Fe respectively then their blueshift suggests that they are the 
ionised signature of the high velocity outflowing wind previously 
detected by Pounds et al. (2004a) and Terashima et al. (2009). We 
are able to model the absorption lines with a single absorber with a 
turbulence velocity of a = 3 000 km s -1 and best-fitting parame- 
ters of log ^ ~ 4.1 and Nn ~ 8.4 x 10 22 cm~ 2 . The zone appears 
to correspond to an outflow velocity of v ont ~ — 5 800 km s -1 
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(~— 0.02c) and is perhaps the highly ionised signature of the high 
velocity outflowing zones of gas detected by Collinge et al. (2001) 
and Steenbrugge et al. (2009). We constrain the location of the ma- 
terial to be on the order of a few 1-d from the black hole, well 
within the dusty torus for this object and may originate in an ac- 
cretion disc wind perhaps depositing up to 10 56 -10 57 erg through- 
out the lifetime of the AGN. We note that we do not, however, 
find any requirement in these data for the higher velocity outflow 
(vout ~ — 30 000 km s _1 ) at Fe K reported by Pounds & Vaughan 
(2010) in their results of a recent XMM-Newton RGS observation. 

Finally, the long-term spectral variability can be modelled 
simply by allowing for changes in the normalisation of the intrinsic, 
unabsorbed power-law component with respect to a quasi-constant 
hard reflection component. This can be interpreted as being due to 
changes in the covering fraction of a Compton-thick absorber ob- 
scuring the intrinsic continuum emission. This is consistent with 
the findings of Miller et al. (2010) who, through studying the ef- 
fects of reverberation in the hard band, find that a global covering 
factor of > 40 per cent of reflecting material is required in this 
source. The constancy of the flux of the near-neutral Fe Ka com- 
ponent across epochs suggests that it does not respond to changes 
in the overall spectral shape and this is also perhaps supportive of 
a partial-covering scenario as opposed to intrinsic changes in the 
continuum flux. 
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Table Al. Table showing the best-fitting parameters of the absorption and emission lines which we detect in the HETG spectrum <2keV. We note that several 



of the absorption 


lines appear to be spectrally resolved with FWHM corresponding to a few hundred to 


~ 1 000 km s 1 . AC corresponds to the change in the 


C-statistic upon modelling the feature in the data. All parameters are given in the rest frame of the host 


galaxy. See Sections 4.1 and 4.2 for further details. 
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Table A2. Table showing the spectral parameters in the rest frame of the best-fitting model simultaneously fit to the 2005 (obs. 1) and 2008 (obs. 2 and obs. 3) 
broad-band 0.6-50.0 keV Suzaku data described in Section 6.5. "Primary power-law continuum: T, photon index; Normalisation. 6 Absorbed power-law. 
c Partial-coverer: Ah, column density; ?, ionisation parameter; / CO v, covering fraction expressed as a percentage. d Neutral reflector: Ap e , iron abundance 
with respect to Solar. e Zone of absorption. The redshift of the zone was fixed at the same value as that of the host galaxy (i.e. z = 0.002336). ^Ionised 
reflector. 9 Highly ionised absorption: Vout, velocity shift. h FeKa emission: Ey lnc , centroid energy; a, intrinsic width; -Fn nc , line flux; EW, equivalent width. 
The symbol p signifies that the parameter has pegged at the maximum / minimum value allowed by the model and the symbol I signifies that the parameter 
is linked to all parameters with the corresponding number. Note that the Fe Ka emission line parameters were obtained by replacing the neutral REFLIONX 
component with a PEXRAV component and a Gaussian. 
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FeKa h 


Eline (keV) 
<r(eV) 
^lino (phcm- 2 s- 
EW(eV) 




6.40i° ° 2 

71 + 19 
11 -23 

2.06±oi 2 
201^22 


6 40+ ' 02 
D -^ u -o.oi 

68 ±33 
2 46+ ' 41 

^•^ u -0.36 


6.41 ±0.02 
< 55.3 
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